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ABSTRACT 


Methods  and  procedures  have  been  developed  for  extrapolating  ground-test 
aerothermodynamic  data  to  flight  conditions  for  four  geometric  elements.  These 
elements  are  a  sharp  unyawed  delta  wing,  deflected  control  surfaces  in  com¬ 
pression  and  expansion,  surface  roughness  elements  typified  by  shallow  convex 
surface  waves  and  rectangular  grooves,  and  a  yawed  circular  cylinder  inter¬ 
secting  a  sharp  flat  plate.  Graphs  and  charts  are  provided  for  rapid  numerical 
computation  of  factors  to  be  used  in  extrapolating  the  wind-tunnel  data  for  these 
geometric  elements  to  flight  conditions.  An  evaluation  of  the  simulation  capa¬ 
bilities  of  operational  hypersonic  ground  facilities  is  also  included. 
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SECTION  I 


INTRODUCTION 


Full-scale  free-flight  tests  of  maneuverable  reentry  vehicles  and  hypersonic 
cruise  vehicles  are  highly  desirable,  since  this  is  the  only  way  in  which  all  of  the 
Row  conditions  influencing  design  parameters  can  be  encountered  simultaneously. 
In  most  cases,  the  cost  of  such  an  approach  is  prohibitive.  As  a  result,  the  ma¬ 
jority  ot  all  hypersonic  aerothermodynamic  data  used  in  configuration  design  are 
obtained  in  ground  facilities  (i.e. ,  wind  tunnels,  shock  tunnels,  hotshot  tunnels, 
helium  tunnels,  and  shock  tubes)  that  cannot  simultaneously  simulate  all  of  the 
conditions  encountered  in  flight.  The  effects  of  simulation  deficiencies  in  ground 
lacilities  on  the  observed  aerothermodynamic  characteristics  of  maneuverable 
reentry  configurations  are  discussed  in  references  1  and  2.  Hypersonic  cruise 
vehicles,  though  not  subjected  to  as  severe  an  environment  as  the  reentry  vehicle 
will  encounter  some  partial  simulation  problems,  particularly  in  regions  of  mu¬ 
tual  flow- field  interference  (e.g. ,  at  the  engine  cowl  lip,  or  at  wing-body  control 
surface  intersections). 

Successful  design  of  hypersonic  vehicles  is  dependent  upon  the  development 
of  analytic  methods  and  procedures  that  can  account  for  the  lack  of  simulation  in 
ground- test  facilities  and  enable  confident  application  of  ground-test  data  to  re- 
entry  vehicle  design.  This  report  presents  an  evaluation  of  the  simulation  capa¬ 
bility  of  ground  facilities  available  for  development  testing  of  scaled  reentry  and 
ypersonic  cruise  vehicles.  Aerodynamic  parameters  such  as  Reynolds  number 
Mach  number,  Lewis  number,  and  Knudsen  number  have  been  evaluated  for 
ground-facility  stream-flow  conditions  and  are  presented  in  Section  II.  These 
parameters  are  used  in  conjunction  with  velocity,  total  enthalpy,  and  free-stream 
temperature  to  illustrate  the  simulation  limitations  of  these  facilities.  In  the 
present  analysis,  the  characteristics  and  the  simulation  capability  of  hyper¬ 
velocity  gun  facilities  have  not  been  included  since  very  few  heat-transfer  and 
pressure  distribution  investigations  using  complex  configurations  have  been  ob¬ 
tained  in  this  type  of  facility. 

Extrapolation  methods  and  orocedures  for  the  aerothermodynamic  analysis 
of  four  geometric  elements  are  presented  in  Section  III.  These  elements  include 
a  8  arp  unyawed  delta  wing,  deflected  control  surfaces  in  compression  and  ex¬ 
pansion,  surface  roughness  elements  typified  by  shallow  convex  waves  and  rec- 
tangu  ar  grooves,  and  a  yawed  circular  cylinder  intersecting  a  sharp  flat  plate. 
Extrapolation  charts  and  curves  are  presented  in  this  section  that  give  numerical 
results.  Extrapolation  factors  are  presented  directly,  or  a  combination  of  fac- 
ors  is  presented  that  will  allow  simple  formulation  of  the  extrapolation  factor. 

The  use  of  the  charts,  and  of  the  extrapolation  procedure  is  also  explained  in 
this  section.  The  flight  conditions  are  represented  by  a  reentry  corridor  defined 
y  a  wing  loading  parameter,  W/SCL,  range  between  50  and  500  and  by  a  flight 
envelope  between  10'2  and  10'6  atmospheres  at  a  free-stream  Mach  number  of  10. 
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The  theoretical  method  used  to  determine  sharp-delta-wing  centerline  heat- 
transfer  is  outlined  in  Appendix  A.  A  new  flow- field  method  was  developed  to 
determine  real-gas  and  leading-edge  sweep-angle  effects  on  flow-field  gas  prop 
erties  at  the  centerline  of  a  sharp  delta  wing.  This  method  is  presented  in 
Appendix  B.  Existing  sharp- delta- wing  spanwise  flow-field  methods  for  angles 
of  attack  less  than  about  15  degrees  (reference  18)  and  gx-eater  than  about  45 
degrees  (reference  36)  have  been  extended  to  the  intermediate  angle  of  attack 
range  between  these  limits.  The  extended  method  bridges  the  analytic  methods 
shown  in  references  18  and  36,  and  postulates  a  parting  line  theoiy  for  some 
conditions  in  the  intermediate  angle  of  attack  range  that  predicts  maximum  heat 
transfer  will  occur  on  the  delta  wing  surface  inboard  of  the  leading  edge  and 
outboard  of  the  centerline.  This  method  is  pi-esented  in  Appendix  C. 
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SECTION  II 


GROUND- FACILITY  FLIGHT  CONDITION  SIMULATION  CAPABILITY 


In  the  past  few  years,  continual  improvement  in  material  capability  and 
cooling  techniques  has  allowed  wind-tunnel  designers  to  extend  the  thermal 
energy  and  maximum  pressure  capability  of  wind  tunnels,  shock  tunnels,  and 
hotshot  tunnels.  Nevertheless,  total  enthalpy  deficiencies  still  exist  in  these 
types  of  facilities  which  affect  flow- field  and  boundary- layer  characteristics, 
heat  transfer,  and  non-dimensional  parameters  such  as  Reynolds  number, 

Prandtl  number,  and  Lewis  number. 

Hypervelocity  gun  facilities  are  being  used  to  simulate  high  enthalpy  effects 
on  models  in  a  free-llight  environment.  The  models  generally  are  very  small 
and  simple,  however,  and  extensive  local  measurements  of  pressure  and  heat 
transfer  are  not  possible  at  this  time.  The  present  analysis  does  not  include 
this  type  of  ground  facility. 

A  discussion  of  the  frec-stream  environment  that  can  be  obtained  in  air- 
nitrogen  tunnels  having  test-section  diameters  on  the  order  of  four  feet  or  greater 
and  in  helium  tunnels  with  test-section  diameters  of  about  two  feet  is  presented 
in  this  section.  The  operating  conditions  for  air  and  nitrogen  tunnels  considered 
in  the  present  analysis  are  characteristic  of  Tunnels  B,  C,  and  F  at  Arnold  Engi¬ 
neering  Developmental  Center  (reference  3)  and  of  the  shock  tunnel  at  Cornell 
Aeronautical  Laboratory  (reference  4).  Operating  characteristics  for  helium 
tunnels  pertain  to  the  facility  at  NASA- Langley1. 


Shock-tube  diameters  do  not  approach  the  test-section  diameters  of  the  facil¬ 
ities  described  above.  Their  range  of  operation,  however,  is  presented  for  com¬ 
parison  with  the  range  of  operation  in  air,  nitrogen.and  helium  facilities.  The 
environmental  range  of  shock  tubes  has  been  determined  for  the  present  analysis 
using  data  from  reference  5. 


Effects  of  total  enthalpy  deficiencies  on  some  nondimensional  parameters 
are  presented  in  this  section.  A  qualitative  summary  of  real-gas  effects  on  free- 
stream  gas  properties  in  ground- facility  test  sections  is  also  presented. 


1  These  characteristics  were  obtained  during  a  facility  tour  at  the  25th  Super¬ 
sonic  Tunnel  Association  Meeting  at  NASA- Langley  Research  Center,  Langley 
Field,  Virginia,  May  1966. 
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1.  GROUND- TEST  FACILITIES 


a.  Energy  Simulation 

The  ability  to  simultaneously  simulate  all  of  the  parameters  important  to 
hypersonic  reentry  in  ground  facilities,  and  to  provide  reliable  aerothermody- 
namic  data  is  discussed  in  references  1  and  2.  It  is  shown  in  those  two  refer¬ 
ences  that  hypersonic  ground  facilities  which  employ  air  or  nitrogen  as  a  test 
medium  cannot  simulate  the  energy  levels  associated  with  reentry  vehicles. 

Shock  tubes  can  achieve  flow  velocities  equivalent  to  high-speed  flight  conditions 
but  at  the  expense  of  a  loss  in  Mach  number  and  free-stream  gas  property  simu¬ 
lation.  Helium  tunnels,  on  the  other  hand,  have  the  advantage  of  achieving  hyper¬ 
sonic  Mach  numbers  with  stagnation  temperatures  on  the  order  of  ambient  tem¬ 
perature.  The  helium  facility,  however,  simulates  neither  the  high  enthalpies, 
nor  the  real-gas  effects  associated  with  flight. 

These  conditions  are  illustrated  in  figure  1  in  terms  of  free-stream  Mach 
number,  M„,  and  free-stream  velocity,  V„.  In  this  figure,  free-stream  Mach 
number  is  used  as  a  measure  of  compressibility  effects.  Duplication  of  this 
parameter  is  required  to  insure  a  realistic  shock-wave  geometry.  Free-stream 
velocity  is  indicative  of  the  kinetic  energy  of  the  stream.  This  figure  indicates 
clearly  that  present  ground  facilities  cannot  simultaneously  duplicate  the  high 
Mach  number  and  velocity  encountered  in  the  region  of  most  critical  reentry 
heating,  which  occurs  within  the  orbital  reentry  corridor  at  velocities  between 
about  17,000  a r  '  22,000  fps  and  at  Mach  numbers  between  16  and  23. 

The  kinetic  and  thermal  energy  capability  of  ground  facilities  provides  a 
criterion  that  allows  separation  of  these  facilities  into  three  sub-groups.  Hyper¬ 
sonic  air  and  nitrogen  tunnels  form  the  first  group  in  which  a  bulk  of  all  aero- 
thermodynamic  design  testing  is  accomplished.  Hypersonic  helium  tunnels  form 
the  second  group,  and  shock  tubes  the  third.  The  energy  capability  of  these  three 
groups  of  tunnels  with  respect  to  the  energy  available  during  reentry  is  best 
illustrated  by  a  figure  which  shows  the  relationship  between  total  energy  and 
thermal  energy.  This  is  done  in  figure  2  where  total  enthalpy,  iQ,  representing 
the  sum  of  the  kinetic  and  thermal  energy  of  the  flow,  is  shown  in  terms  of  the 
free-stream  temperature,  Ta. 

The  energy  capability  of  hotshot  tunnels,  shock  tunnels,  and  conventional 
wind  tunnels  is  illustrated  by  the  region  in  the  center  of  the  figure.  This  region 
is  bounded  at  low  free-stream  temperatures  by  the  onset  of  gas  liquefaction  in 
the  test  section.  This  boundary  is  approximate,  in  that  it  corresponds  to  free- 
stream  pressures  and  temperatures  characteristic  of  the  facilities  being  ana¬ 
lyzed.  The  upper  boundary  is  defined  by  the  total  enthalpy  available  in  these 
facilities.  Combustion-driven  shock  tunnels  reach  energy  levels  in  excess  of 
4,000  Btu/ibm,  whereas,  hotshot  tunnels  and  "regular"  shock  tunnels  (i.e. , 
those  using  inert  gases  for  the  driver  gas)  are  able  to  reach  energy  levels  on 
the  order  of  2,  800  Btu/lbm.  (Combustion-driven  shock  tunnels  are  separated 
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fiom  regular"  shock  tunnels  because  stream- flow  characteristics  for  these 
facilities  are  not  well  defined  for  all  operating  conditions.)  Conventional  wind 
tunnels  can  obtain  total  enthalpies  on  the  order  of  1,000  Btu/lbm.  It  is  indicated 
in  this  figure  that  the  flight  environment  can  be  simulated  by  the  hotshot  and 
shock  tunnels  up  to  a  free-stream  velocity  of  about  11,000  fps. 

Extremely  low  free-stream  static  temperatures  (i.e. ,  approximately  2°  to 
10°  Rankine)  are  characteristic  of  helium  tunnels.  A  lower  limit  is  imposed  on 
the  free-stream  temperature  that  is  based  on  liquefaction.  Total  enthalpies 
above  ambient  levels  can  be  obtained  in  this  type  facility  by  heating  the  helium. 
However,  the  majority  of  present-day  operational  helium  facilities  operate  with 
total  enthalpies  which  correspond  to  ambient  temperatures. 

Shock  tubes  can  reach  energy  levels  on  the  order  of  30,000  Btu/lbm  and 
velocities  on  the  order  of  24,000  fps.  This  group  of  facilities,  however,  does 
not  lend  itself  to  development  tests  designed  to  obtain  complex  shape  aerodynamic 
heat-transfer  and  pressure  measurements. 

b.  Reynolds  Number  Simulation 

The  capability  to  simulate  the  free-stream  Reynolds  number  in  hypersonic 
air,  nitrogen,  and  helium  tunnels  is  shown  in  figure  3  in  terms  of  free-stream 
Mach  number.  The  free-stream  Reynolds  number  is  based  on  a  one- foot  model 
in  air  and  nitrogen  tunnels,  and  a  0.  25  foot  model  in  helium  tunnels.  The  region 
of  free-stream  Reynolds  numbers  associated  with  reentry  is  also  shown.  This 
region  is  bounded  by  values  corresponding  to  trajectories  with  a  wing  loading 
parameter,  W/SCL,  range  between  50  and  500.  The  region  of  free-stream 
Reynolds  numbers  pertaining  to  sustained  hypersonic  cruise  between  altitudes 
of  100,  000  and  300,  000  feet  at  M„  =  10  is  also  shown.  The  free-stream  Reynolds 
numbers  for  flight  were  computed  assuming  a  prototype  30-foot  vehicle.  Simu¬ 
lation  oi  free-stream  Reynolds  number  and  Mach  number  over  these  two  flight 
regimes  is  shown  to  be  quite  good. 


Simulation  of  the  free-stream  Reynolds  number  and  the  Mach  number  in 
ground  facilities  does  not  insure  that  boundary -layer  development  over  the  test 
model  will  be  similar  to  boundary- layer  development  over  the  prototype  vehicle 
in  flight.  Boundary- layer  development  is  dependent  upon  the  boundary- layer 
edge  Reynolds  number,  Ree,  the  edge  Mach  number,  Me,  and  the  ratio  of  the 
wall  temperature  to  the  boundary- layer  edge  temperature,  Tw/Te.  Large  dif¬ 
ferences  between  wind-tunnel  and  flight  values  for  Re  and  M„  have  been  calcu¬ 
lated  for  sharp  flat  plates.  Differences  in  Ree  and  Me  are  illustrated  for  two 
free-stream  Mach  numbers  in  figure  4  in  terms  of  angle  oi  attack,  Q. 

Edge  Reynolds  number  differences  depend  also  somewhat  on  the  geometry 
as  is  illustrated  by  figure  5  for  sharp  flat  plates  and  sharp  cones.  In  this  figure, 
the  Reynolds  number  ratio  of  figure  4  has  been  inverted  and  presented  as  a  func¬ 
tion  of  a  and  the  ratio  of  total  enthalpies  iWT/iFLT.  The  band  indicates  the 
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effect  of  free -st  re  am  pressure  in  the  range  from  10_1  to  10-6 
at  M.  =10  and  from  10"®  to  10-®  atmospheres  at  M,,  =  20. 
enthalpy  ratio,  the  differences  in  Ree  at  M,  =  20  approach  a 
high  angles  of  attack  for  both  geometries. 

c.  Lewis  Number  and  Knudsen  Number  Simulation 

1)  Lewis  Number 

The  dimensionless  group  (P<L,f  D/k^ )  Is  defined  as  the  Lewis  number,  Le. 
It  represents  the  ratio  of  energy  transported  by  diffusion  to  that  transported  by 
conduction.  The  subscript  f  on  the  specific  heat,  Cp,  and  thermal  conductivity, 
k,  indicate  that  these  parameters  represent  a  frozen  gas  composition  (i.e. ,  no 
chemical  reactions  are  taking  place).  The  other  two  terms  in  this  group  are  the 
density,  P ,  and  D,  the  binary  diffusion  coefficient  between  the  dissociated  and 
undissociated  species.  This  parameter  is  generally  used  as  a  correction  factor 
that  accounts  for  diffusion  effects  on  convective  heat  transfer.  For  instance,  in 
reference  6,  it  appears  in  the  correlation  equation  for  stagnation-point  heating 
as 


atmospheres 
For  the  lowest 
factor  of  10  at 
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whereas,  in  reference  7,  stagnation -point  heat  transfer  for  Lewis  numbers 
different  from  unity  is  represented  by  the  expression 
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In  equations  (1)  and  (2),  ijj  is  the  enthalpy  absorbed  in  the  process  of  dissociation, 
iQ  is  the  total  enthalpy  of  the  undissociated  gas,  and  w  refers  to  conditions  at  the 
wall.  It  is  shown  in  reference  7  that  for  a  highly  cooled  wall  (i.e. ,  Tw/T  «  0. 1), 
equations  (1)  and  (2)  agree  to  within  a  few  percent. 


Lewis  number  comparisons  can  be  made  between  flight  and  ground  facilities 
that  use  air  or  other  multicomponent  gases  as  a  test  medium.  The  range  of 
Lewis  numbers  that  can  be  obtained  in  facilities  using  air  as  a  test  medium  is 
presented  in  Table  I.  Curves  of  reference  8  indicate  that  the  Lewis  number  can 
range  from  about  0.6  to  1.5  in  flight  and  will  generally  be  less  than  unity  for 
that  portion  of  the  reentry  corridor  where  heating  is  most  severe. 


TABLE  I 

GROUND  FACILITIES 

Wind  tunnels,  shock  tunnels, 
and  hotshot  tunnels 
Shock  tubes 
FLIGHT 


Lewis  Number  Range 
1.0  -  1.5 

0.6  -  1.5 
0.6  -  1.5 
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Values  in  Table  I  indicate  that  Lewis  number  for  the  entire  range  of  reentry 
can  only  be  simulated  in  shock  tubes. 

Since  the  Lewis  number  is  only  used  for  small  corrections^  to  the  heat- 
transfer  calculation,  its  duplication  is  not  as  important  as  is  the  duplication 
of  the  other  parameters  already  discussed. 

2)  Knudsen  Number 

The  Knudsen  number,  Kn,  is  defined  by  the  dimensionless  parameter  X/f  , 
where  A  is  the  mean-free-path  of  the  moleculeB  comprising  the  free-stream 
fluid,  and  f  is  a  characteristic  body  length.  The  Knudsen  number  is  generally 
used  as  a  criterion  to  determine  if  a  fluid  can  be  treated  as  a  continuous  medium. 
WhenKn«l,  it  is  assumed  that  a  continuum  behavior  exists.  Free-stream 
conditions  in  the  three  groups  of  facilities  considered  in  this  anulysis  meet  this 
criterion. 


2.  REAL-GAS  EFFECTS  IN  GROUND-FACILITY  NOZZLE  FLOW 

Real -gas  effects  on  ground -facility  nozzle  flow -field  properties  are  directly 
related  to  chemical  changes  that  result  from  high  stagnation  chamber  tempera¬ 
tures.  For  instance,  up  to  temperatures  on  the  order  of  1200  “R  air  can  be  con¬ 
sidered  to  behave  as  a  perfect  gas.  At  temperatures  greater  than  about  1200  'R, 
appreciable  vibration  of  the  gas  molecules  begins,  but  no  chemical  reactions 
occur.  The  region  of  vibrations  extends  from  about  1200  ’R  to  about  3600  "R, 
Dissociation  of  the  oxygen  contained  in  air  molecules  becomes  appreciable  at 
about  3600  °R.  Dissociation  of  the  nitrogen  contained  in  air  1b  the  next  major 
chemical  reaction  that  occurs.  At  about  10,  000  “R,  ionization  of  the  atoms  be¬ 
comes  significant,  and  large  quantities  of  positively  and  negatively  charged 
particles  appear  in  the  air.  As  this  temperature  is  increased,  the  fluid  becomes 
more  and  more  like  a  plasma  with  an  electrical  conductivity  comparable  to  sea 
water.  Nitrogen  and  helium  react  to  high  temperatures  in  a  similar  manner. 

Total  enthalpy  levels  in  present-day  ground  facilities  are  high  enough  to 
cause  significant  real-gas  effects  on  gas  properties.  Real-gas  effects  on  helium 
tunnel  nozzle  flows  are  generally  small  due  to  the  low  energy  state  in  the  stag¬ 
nation  chamber.  Air-nitrogen  tunnels,  on  the  other  hand,  operate  at  stagnation 
temperatures  high  enough  to  cause  molecular  dissociation  in  the  stagnation  cham¬ 
bers.  For  example,  it  is  shown  in  reference  9  that  for  air  in  chemical  equilib¬ 
rium  at  a  pressure  of  one  atmosphere  and  a  temperature  of  5900  °R  the  oxygen 
contained  in  the  air  is  approximately  50%  dissociated.  As  a  result  of  this  dis¬ 
sociation,  real -gas  effects  on  test -section  gas  properties  of  the  hypersonic 


2  The  maximum  effect  of  diffusion  on  heat  transfer  with  Le  =  1.4  in  an  air 
facility  is  about  5%.  In  flight,  at  velocities  in  excess  of  18,  000  fps  the  effect 
of  Le  is  approximately  10%  at  the  stagnation  point. 
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free  stream  are  dependent  upon  not  only  enthalpy  and  entropy,  but  also  upon  the 
time  required  for  atom  recombination  with  respect  to  the  time  required  for  ex¬ 
pansion,  The  degree  of  dissociation  must  be  nearly  zero  in  the  free  stream  at 
the  test  section  in  order  to  simulate  the  free  stream  that  will  be  encountered  at 
most  flight  conditions. 

Since  recombination  of  a  dissociated  gas  requires  a  finite  amount  of  time, 
it  is  possible  that  the  stream  flow  may  reach  the  test  section  partially  dissociated. 
Expanding  flows  in  chemical  nonequilibrium  are  difficult  to  assess  analytically 
but  can  be  evaluated  for  the  limiting  cases  of  equilibrium  and  frozen  expansions 
(e.g. ,  references  1,  10,  and  11).  If  the  flow  is  considered  frozen  as  it  expands, 
it  is  most  probable  that  the  flow  will  tend  toward  equilibrium  again  as  a  result 
of  compression  near  a  test  model.  This  is  the  inverse  of  what  happens  in  flight 
where  dissociation,  if  it  occurs,  takes  place  behind  the  shock  wave  generated  by 
a  hypersonic  vehicle  or  in  the  boundary-layer  near  the  surface.  As  a  result, 
stream-flow  gas  properties  typical  of  flight  conditions  can  only  be  simulated  in 
air-nitrogen  tunnels  when  the  stream  is  in  equilibrium.  It  is  shown  in  reference 
10  that  for  stagnation-chamber  pressures  greater  than  about  35  atmospheres, 
and  total  enthalpies  less  than  about  2650  Btu/lbm,  the  expansion  of  air  in  a  nozzle 
is  essentially  in  equilibrium. 

Charts  are  provided  in  reference  1  for  determining  real -gas  effects  in  hyper¬ 
sonic  nozzle  flows  in  equilibrium  air  from  knowledge  of  stagnation-chamber  con¬ 
ditions.  A  more  detailed  study  of  these  effects  has  been  presented  in  reference 
12  with  results  that  allow  rapid  determination  of  test -section  gas  properties. 
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SECTION  in 


EXTRAPOLATION  METHODS 


This  section  describes  methods  and  procedures  for  extrapolating  wind-tunnel 
heat -transfer  data  to  flight  conditions.  The  term''extrapolation”  denotes  adjust¬ 
ment  of  wind-tunnel  data  to  account  for  lack  of  complete  simulation  of  flight  con¬ 
ditions  in  the  air-nitrogen  and  helium  tunnels  discussed  in  Section  n. 

A  series  of  charts  is  presented  that  can  be  used  in  the  development  of  extrap¬ 
olation  factors  for  the  following  geometric  elements;  (1)  sharp  delta  wings, 

(2)  deflected  control  surfaces,  (3)  interference  geometries  typified  by  yawed 
cylinders  attached  to  a  sharp  flat  plate,  (4)  shallow  convex  surface  waves  and 
rectangular  grooves.  These  charts  are  evaluated  at  conditions  representative  of 
both  the  flight  environment  defined  by  the  reentry  corridor  illustrated  in  figure  6 
and  the  environment  in  air-nitrogen  and  helium  tunnels.  (Some  operational  air- 
nitrogen  tunnels  can  simulate  the  environment  of  sustained  hypersonic  cruise. 

As  a  result,  extrapolation  factors  are  not  generally  required  for  heat -transfer 
and  pressure  measurements  made  in  this  flight  regime  and  have  not  been  specif¬ 
ically  treated  in  this  report. ) 

Within  the  reentry  corridor  shown  in  figure  6,  dependence  of  extrapolation 
factors  on  pressure  is  on  the  order  of  3%;  hence,  average  pressures  are  used  to 
represent  the  flight  regime.  Changes  in  the  extrapolation  factors  that  are  due 
to  variations  in  total  enthalpy  at  the  same  Mach  number  in  flight  are  on  the  same 
order  and  have  been  neglected  in  this  report.  All  charts  representing  flight  con¬ 
ditions  are  calculated  for  a  free-stream  temperature  equal  to  400  °R,  and  a  wall 
temperature,  Tw,  equal  to  2000  °R,  Wind-tunnel  charts  are  calculated  for 
Tw  540  "R,  All  calculations  are  made  for  flows  in  chemical  equilibrium. 
Nonequilibrium  effects  on  the  extrapolation  factors  are  discussed  in  Section 
m-6.  All  calculations  are  made  using  the  thermodynamic  gas  properties  of 
reference  13. 

The  basic  approach  of  the  present  analysis  in  extrapolation  factor  develop¬ 
ment  is  to  use  analytic  methods.  In  this  manner,  those  parameters  influencing 
the  extrapolation  factors  can  be  evaluated  systematically. 

1.  BASIC  EXTRAPOLATION  PROCEDURE 

Application  of  ground-test  data  to  flight  conditions  is  accomplished  in  four 
parts.  First,  the  data  are  extrapolated  from  wind-tunnel  to  flight  enthalpy  con¬ 
ditions  at  constant  Mach  number.  The  differences  between  these  two  conditions 
are  expressed  by  a  factor  Fp  the  subscript  i  denoting  the  fact  that  enthalpy  is 
the  varying  quantity.  In  functional  form,  Fj  appears  as: 
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(3) 


Fi  (hFLT/  Vt'm.,,  WT 

where  h  is  a  normalized  laminar  heat -transfer  coefficient  and  the  subscripts  WT 
and  FLT  denote  wind-tunnel  and  flight  conditions  at  the  wind  tunnel  Mach  number, 
Mw.  WT.  Values  for  Fj  are  shown  in  charts  for  each  of  the  elements. 


The  second  part  of  the  extrapolation  procedure  consists  of  "adjusting"  the 
partially  extrapolated  data  to  account  for  differences  between  the  wind-tunnel 
free-stream  Mach  number  and  the  flight  value.  These  differences  are  expressed 
in  terms  of  the  factor  Fj^,  where  the  subscript  M  denotes  that  the  Mach  number 
is  the  varying  quantity.  The  factor  F^j  is  computed  by  foiming  a  ratio  of  the 
normalized  heat -transfer  coefficients  at  the  flight  Mach  number  to  that  corres¬ 
ponding  to  the  wind-tunnel  Mach  number  as  follows: 

FM  =  (^M^,  FLT^M,,,,  WT  )  (4) 


A  third  extrapolation  factor  Fx  is  required  to  account  for  length  differences 
on  the  sharp  delta  wing.  For  laminar  flow  Fx  is  defined  as 

1/2 

Fx  =  **WT/x  FLT* 

where  x  is  the  distance  measured  from  the  apex  or  leading  edge  of  the  wing. 


A  fourth  extrapolation  factor  which  only  appears  in  the  sharp  delta-wing 
extrapolation  procedure  accounts  for  cross -flow  pressure  gradients  and  is  de¬ 
fined  as 


=  0*  FLT 
=  0*  WT 


(6a) 


This  procedure  is  illustrated  in  the  following  sketch  for  laminar  boundary- 
layer  flow 
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For  the  sharp  delta  wing  in  turbulent  flow,  a  fifth  extrapolation  factor 
Fgf,  is  required.  It  is  based  on  differences  in  the  local  Reynolds  number  and 
is  defined  as 


hL*  FLT 
W  WT 


where  h  ,h 

T/  L 


is  defined  by  equation  (A  19)  in  Appendix  A. 


(6b) 
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2.  SHARP  DELTA  WINGS 


a.  Laminar  Flow 


To  use  the  procedure  illustrated  on  the  previous  page,  it  is  necessary  to 
normalize  wind-tunnel  heat -transfer  data  measured  locally  by  a  laminar  refer¬ 
ence  heat -transfer  coefficient,  href,  defined  as  the  stagnation-point  heat -transfer 
coefficient  for  a  hemisphere  with  a  one-foot  radius.  Wind-tunnel  stagnation- 
point  heat -transfer  data  can  be  transformed  to  hrej  by  the  expression 

1/2 

^refHvT  "  (ho)WT<r)  (7) 

where  r  is  the  radius  of  the  test  model  in  feet.  Local  heating-rate  data  is 
normalized  by  href  for  the  free-stream  Mach  number  at  which  the  data  were 
measured  to  form  the  ratio  (  h  /  h  _) 

L  ref  M  WT 

The  extrapolation  factors  at  the  centerline  of  the  delta  wing  are  defined  by 
the  following  expressions 


(8) 

(9) 

(10) 


where  the  subscript  (M<=,  x)WT  denotes  evaluation  at  the  Mach  number  and 
distance  in  the  wind  tunnel  and  the  subscript  (i,x)FLT  denotes  evaluation  at  the 
enthalpy  and  distance  in  flight. 

In  equations  (8)  and  (9)  the  ratio  (h/href)  is  determined  at  the  centerline  using 
the  heat -transfer  method  presented  in  Appendix  A.  Flow-field  gas  properties  re¬ 
quired  to  compute  the  heating  rates  in  proximity  to  the  centerline  are  computed 
using  the  method  described  in  Appendix  B. 


For  the  delta  wing,  the  extrapolation  factor  Fj  in  air-nitrogen  wind  tunnels 
has  been  evaluated  in  terms  of  the  enthalpy  ratio  (ioXVT/^olFLT  anc*  free-stream 
Mach  number,  M^.  It  is  presented  in  figure  7  for  delta  wings  with  leading-edge 
sweep  angles  between  70  and  80-degrees  in  the  angle -of -attack  range  between  10 
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and  60-degrees.  This  factor  for  helium  tunnels  is  presented  in  figure  8  in  terms 
of  the  free-stream  Mach  number  and  angle  of  attack. 

Solutions  for  the  heat -transfer  ratios  required  to  form  the  extrapolation 
factor  Fjvi  are  presented  in  figure  9.  The  ratio  (h/href)  for  flight  is  obtained 
at  flight  Mach  numbers.  Figure  9  is  presented  for  delta  wings  having  leading- 
edge  sweep  angles  of  70,  75,  and  80-degrees. 

The  extrapolation  factor  Fx  is  obtained  using  equation  (10). 
b.  Turbulent  Flow 


Extrapolation  of  turbulent  heat -transfer  data  to  flight  conditions  is  accom¬ 
plished  using  an  extension  of  the  laminar  method  described  above.  This  proced¬ 
ure  is  diagrammed  in  the  lower  portion  of  the  preceding  sketch.  To  begin  the 
extrapolation  procedure,  the  wind-tunnel  turbulent  heat  transfer  data,  hf,  is 
also  normalized  by  href,  the  laminar  reference  heat -transfer  coefficient  discus¬ 
sed  in  the  preceding  text.  In  addition  to  the  factors  Fj  and  F^,  an  extrapolation 
factor  that  incorporates  differences  in  the  local  Reynolds  number  is  required. 
This  factor,  Fr  ,  is  a  ratio  of  turbulent  to  laminar  heat -transfer  coefficient 
ratios  at  flight  and  wind-tunnel  conditions  and  is  defined  as 

VV.FLT 
Rr~  VhL\vT 


The  ratio  hj/hL  is  evaluated  in  terms  of  the  local  Reynolds  number  by  the 
expression 


h 


0. 185  t/IT 


L  |o.332^1og10  (Rr+3000)J 

nolds  number  c 

p  u  u  x  , 
^r'r  e  eg,  L 


2.584  | 


where  R^  is  the  local  reference  Reynolds  number  defined  as 


R 


2  2 
“a  Fx 


(12) 


(13) 


These  expressions  are  unique  to  the  ”prHr”  heat -transfer  method  described  in 
Appendix  A. 


In  order  to  evaluate  fRr>  the  local  reference  Reynolds  numbers  in  flight  and 
in  the  wind  tunnel  must  first  be  known.  Values  for  Rr  at  the  centerline  of  the 
delta-wing  are  presented  for  wind  tunnel  and  flight  conditions  in  figure  10.  In 
this  figure,  Rr  is  presented  for  x  =  1  foot  and  has  been  normalized  by  (Rr)ref> 
the  Reynolds  number  for  a  sharp  flat  plate  at  a  =  20°  and  x  =  1  foot.  Wind- 
tunnel  values  for  (Rr)ref  are  presented  as  a  function  of  free-stream  Mach  number 
and  temperature  in  figure  11  and  are  normalized  by  the  free-stream  pressure. 
Flight  values  are  presented  in  terms  of  altitude  and  free-stream  velocity  in 
figure  12.  Local  values  for  Rr  are  determined  using  the  expressions 
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[^rejpLT.x  =  1  ft  '  ^ref.FLT  ‘  *FL1 

^rW  [(Rr)refJWT  x  =  lft  •  (Rr)ref,VVT  *  Vt 


(14) 


(15) 


ofthJteflfrTlatI°niaCt0r  FRr  18  Sh°Wn  in  flgure  13  as  a  function  of  a  ratio 

treVcr^r^  ynoids  number  at  mght  conditi°ns  i° that  at  wind- 

The  complete  extrapolation  procedure  that  is  used  to  extrapolate  sham  delta 
wing  centerline  heat -transfer  data  In  turbulent  flow  is  illustrated  below: 


Normalized 
Wind-Tunnel  Data 


0 


X 


c*  Off -Centerline 


iTipfliii 

fv  - 1 1.F4 


The  extrapolation  procedure  for  laminar  and  turbulent  heat -transfer  data 
measured  at  off-centerline  (spanwise)  locations  on  a  delta  wing  consists  of  multi 
Plying  the  extrapolated  centerline  data  by  an  additional  factor,  F^,  defined  as 


(h  /h  .  \ 

<t>  0  =  Q'FT.T 


*  tu/h, 


Vn<t> 


0*WT 


(16) 


^!rhe  hf  *S  the,  local  «Panwise  heat -transfer  coefficient  and  h  .  is  the  center- 
"kich  °  C‘“  m“SUred  at  *  the  (“lowing 


Spanwise  heat -transfer  coefficients  are  calculated  for  hoth  wind  *  , 

flight  conditions  by  the  methods  described  in  Appendix  A,  Flow-field  ^^0- 

in  Appendriex  cetermined  referenCe  13  and  the  Aow-field  method  described 
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.•JET?  rxts.'sr ts.  i£rTm  M,ch  nnmber  20  ■» 
jt%  t™  « *“* !°- 

"°‘wt  ^o)  ,  equal  to  0.25.  ^ 

F  LT 

d-  Dissimilar  Flow-Field  FfWt« 

.he  eh«s‘ d°“r,«  ‘„"r.  «ir‘"s  ""’d-,un"ei  <"*•  *»  •*> «h 

rounding  ,he  de"' wfng  “  TL?'*  '?***“  Wt*”  th«  "»»  «•>-  ~r- 

ofdl.uimllsr  now  fluid,  on  .harp  deitn  wingVort^Z* L^Kudt  “Zj! 

wide  ™£”f  "  giu.'S’S  •“,?“■  *  ,h,rP  delta  «•<"  “  huper.onic  for  . 
field  i.  .ubeonio  .  n  ,  H0W,er'  “  very  high  angle,  of  .flack  the  flow 

field  can  be  tran.mit.ed  ^ream*  and'in  .hi’™'*”0,"'  <i“,urb*noe  ln  the  flow 
flow  field  Snrh  n  Hi  t  k  1  n  tb*s  way  a^ter  the  entire  lower  surface 

vertica!  fins  etc  may  be  caused  b>’  defle^ed  control  surfaces 

IriTy  due  to  ;  £.  of  ,^1  flow  ln  heat  transf-  «*«  5i»- 

Of  course  the  flow  field  Hiff  eld  simulation  rather  than  enthalpy  simulation. 

enthalf^^but'the'res^idtant^effects168  T?  ^  lnduCed  by  dlfferenc -  *-  total 
cannot  be  determined due  to  dis8lmilar  flow  fields 
viscid  subsonic  flow  field  This  ^  develoPed  for  calculating  the  in¬ 

calculating  the  flow  over  a  triwaSlw  w  fhS  approximately  equivalent  to  that  of 

siderably  more  complex  th^  »n^  *  ^  V  at  a"Kle  of  attack-  and  «*°  is  con- 

attack  at  which  the  flow  field  “comes  accompllshed.  The  angle  of 

This  effect  ,s  illustretld  in  fi^reTs  TsZ  real-gaS  effeCt8‘ 

surface  flow  field  on  a  70-deeree  «w»nt  h  B  figure,  it  is  seen  that  the  lower 

20  wind  tunnel  at  a  G0d  a^e *7*3 S  *“7'?*  18  8Ub8°nlC  ln  a  Mach 
same  angle  of  attack.  supersonic  for  a  flight  case  at  the 

3.  DEFLECTED  CONTROL  SURFACES 

age  men  t  ofaer^ynamlc  fTceT  r?''  VehIcles  18  throu«h 

add  flaps.  However  when  7™  3y  t0  accomplish  this  control  is  to 

■he  flap  1.  dl”rS  fl°W  *» 

are  non-existent.  As  a  res^t  Ter^  f  adequately  describe  this  flow  field 
butions  on  deflected  flaps  aerodynamlc  pressure  and  heat -transfer  distri- 
come  from  wind-tunnel  tests.  re8ence  o  upstream  disturbances,  will  largely 

d,.f.t^ »t!Zdtrsr?wrr'  wm-,m  d»*  •»  ~«- 

pr«mded  by  uu,,or„“ZX^ 

these  control  surfaces  are  rienonHu..  186  stream~flow  gas  properties  on 

dynamic  g„  S  e  pr°P'rtlb!  *nd  .hermo- 

,ucface.,  ,he  dynamic  p^  ”“ani  d<,n“,'d  comprea.lon 

y  mic  pressure  and  local  pressure  at  the  hinge  line  exert  the 
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L  I  .  I1161106  °"  extrapolation  factors  presented  in  this  report.  On  ex- 
Sarv  I  ni  rapolatlon  methods  developed  here  are  dependent  upon  the 
h7„^yHr  dl8Placement-thlckness.  pressure,  and  heat -transfer  coefficient 
at  the  hinge  line.  These  extrapolation  methods  are  described  below. 


a.  Compression  Surfaces 


It  isposlib  e  that  th  f  ,nt°  tHe  Stream  Wil1  comPress  the  oncoming  flow. 
Le  aZl  if  S  m  ^  layer  °f  thlS  8treanl  wil1  ^parate  from  the  sur- 

da?a  of  i^feiice  iTwi116'  A"  6Xample  of  thls  Phenomenon  is  shown  by  the 
Sr  '  Whea  ^Paration  occurs,  the  pressure  and  aerodynamic 
nesting  distributions  are  altered  significantly. 

usings  °f  h6at  transfer  throi«h  a  'amlnar  separated  region 

using  an  extension  of  boundary-layer  theory  is  presented  in  reference  15  Thit 

analysis  indicates  that  average  aerodynamic  heat  transfer  is  reduced  50%  as 

16 m  Vd  '  T,  -  be“  V,rl'led  b>'  *“*■  <»“  reference,  ' 
att ihml!;  a  g  regl°ns  °f  flo>V  separation,  their  location  of  re- 

been  pTsslwe  aCC°mpanying  heatin^  distributions  analytically  has  not  yet 

mn^n°ther  approach*  the  one  “sed  in  the  present  analysis,  is  to  predict  the 
more  severe  situation  —  one  resulting  in  maximum  heating  rates.  An  approxi- 
mate  method  for  predicting  an  upper  limit  on  the  local  heat  transfer  to  flips 

bas^i  tK  aminar  °r  tUrbUlent  8treamS  18  Presented  refeSnce  14  "it  is 
based  on  the  expression 


max 


max 


(16a) 


hl  'hl 

where  h  is  the  maximum  heat -transfer  coefficient  expected  on  the  control 
surface  hHL  is  the  heat-transfer  coefficient  at  the  hinge  line,  p"  is  the  pres 
sure  predicted  using  the  oblique-shock  method,  and  PHL  is  the  pSure  at  the 
hinge  line.  This  method,  which  is  verified  by  data  contained  in  that  report  and 

n  references  18  and  19,  forms  the  basis  for  the  extrapolation  procedure  illus¬ 
trated  in  the  following  sketch:  procedure  Ulus- 
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to  both  laminar  and  turbulent  flows 
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To  extrapolate  the  maximum  heating  rate  measured  on  a  compression  surface, 
the  deflected  flap  heat-transfer  data  must  first  be  normalized  using  heating  rates 
measured  just  upstream  of  the  hinge  line.  (It  is  important  to  note  that  the  measured 
heat-transfer  data  at  the  hinge  line  and  on  the  control  surface  must  be  either  lami¬ 
nar  or  turbulent,  not  transitional.)  The  wind-tunnel  data  are  then  extrapolated  to 
flight  conditions  using  the  extrapolation  factor  Fjj^  ,  which  is  defined  as  follows: 

(bmax^HL^l  .FLT 
lM  <hmax/'hHL>M,WT 

where  the  subscripts  M,FLT  and  M,WT  denote  the  flight  Mach  number  and  the 
wind-tunnel  Mach  number,  respectively. 

In  equation  (17),  the  factor  (hmax/hf|L)  is  determined  by  equating  it  to  a  pres¬ 
sure  ratio  determined  using  the  oblique-shock  method,  and  assuming  the  flow  field 
upstream  of  the  hinge  line  is  uns sparated.  The  pressure  ratio  has  beep  correlated 
for  both  wind-tunnel  and  flight  conditions  in  terms  of  the  parameter  (2q/P)jjL, 
where  4  is  the  dynamic  pressure  and  P  is  the  pressure,  both  evaluated  at  the 
hinge  line.  This  correlation  is  presented  in  figure  16  for  air-nitrogen  wind- 
tunnels,  where  the  ratio  P/Phl  has  been  replaced  by  the  ratio  hmax/hHL.  The 
lines  labeled  y  -  1.4  in  figure  16  represent  ideal-gas  solutions  in  air.  For  a 
constant  flap -deflect!  on  angle  6,  all  of  the  solutions  are  confined  to  a  narrow 
region  which  varies  with  (2<i/P)jjl>  and  which  widens  with  increasing  flap- 
deflection  angles.  At  fixed  values  of  6  and  (2^/P)jjl,  the  upper  bound  on  hmax/ 
hpL  is  represented  by  the  ideal-gas  solution.  With  increasing  departure  from 
ideal -gas  conditions,  solutions  shift  toward  the  lower  bound  of  the  region.  The 
lower  bound  in  figure  16 ,  labeled  T  eff  =  1.1,  represents  flight -case  solutions 
at  a  free-stream  Mach  number  of  20  and  a  pressure  of  10'3  atmospheres.  Real- 
gas  effects  in  the  flow  field  fostered  by  the  flap  are  negligible  in  helium  and  fig¬ 
ure  17  presents  heating  rate  ratios  applicable  to  all  helium  tunnels  operating  at 
stagnation  temperatures  equal  to  the  ambient  temperature. 

Values  of  the  parameter  (2q/P)jjL  are  dependent  upon  the  geometry  upstream 
of  the  hinge  line  and  must  be  evaluated  independently.  In  the  correlations  shown  in 
figures  16  and  17,  the  only  restriction  placed  on  the  flow  field  upstream  of  the 
hinge  line  is  that  it  be  two-dimensional  and  uniform. 

b.  Expansion  Surfaces 

Analytic  methods  for  obtaining  viscosity-induced  pressure  and  heat  transfer 
to  a  surface  downstream  of  a  fully  developed  boundary  layer  in  an  expanding  hyper¬ 
sonic  inviscid  flow  field  are  not  well  developed.  This  problem  is  characterized 
in  references  2  and  14  by  the  expansion  of  a  non -rotational  hypersonic  flow  field 
about  a  sharp  corner.  In  these  two  references,  it  was  assumed  that,  (1)  Prandtl's 
boundary-layer  concept  was  valid  such  that  the  viscous  and  inviscid  flow  regions 
could  be  separated,  and  (2)  the  effective  turning  angle  of  the  inviscid  flow  field  was 
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J  ““  gr0W,th  of  viseous  boundary  layer.  Both  of  these  references 

T*?*  Calculatlon  of  *is  il™  field  requires  simultaneous  solu¬ 
tion  of  the  boundary-layer  equations  and  the  equation  for  a  Prandtl-Meyer  expan- 

houn'ri  “  ,  06  2’  an  iteratiVe  Pr0CeSS  based  aP°n  ««!*»  from  a  nonsimuar 

boundary-layer  program  (reference  19)  and  the  Prandtl-Meyer  equation  was  used 

field* was^i'  h n  eH5?re88/°n  f°r  the  6ffeCtlve  angle  Ja  two-dimensiZl Tot' 

field  was  derived  in  reference  14,  using  phyiscal  reasoning  in  conjunction  with  an 

ssumed 1  similarity  for  the  boundary -layer  profiles.  The  effective  turning  angle 
Av,  as  defined  in  that  document,  is  6  g 


in  which 


and 


Au  =  AF  +  ^ 

pii+ 1  dx 


P  =  -2— 
PHL 


(18) 


(19a) 


X  = 


x  -  x, 


HL 


HL 


Ax 

6*HL 


(19b) 


h  (  )  f  ia  the  l0cal  Pressure.  PHL  is  the  hinge  line  pressure, 

f  w  flap  turning  angle,  x  is  the  distance  from  the  leading  edge  to  the  loca- 

is  5. e  disci1  ’  X^k1  dlStanCe  fr°m  th°  Ieadlng  edge  to  line,  6*HL 

is  the  displacement  thickness  at  the  hinge  line,  and  n  is  equal  to  unity. 

Solution  of  equation  (18)  results  in  values  for  the  pressure  ratio  P/Pm  that 

he  forePe?|dtTtf0n>1.y  UPO"  AX/6  HL<  This  suBgests  that  extrapolation  factors  can 
be  formed  that  allow  extrapolation  of  both  pressure  and  heat  transfer  data  for  ex¬ 
pansion  surfaces  to  flight  conditions.  The  procedure  suggested  is  illustrated  below. 
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win“ir„e,a^tl0^Pr0C:dUre;  (h/hHL)WT  and  (P/PHL)WT  are  the  normalized 
wind  tunnel  data.  The  extrapolation  factor  FiM  p  is  defined  as 


F 


<p/phl>flt,  (Ax/6*)ptt 
iM>P  '  <P/PHL)WT,  (Ax/6*)wt 


(20) 


flow  on  the  nln  the  pressure  ratio  (P/PHL)  is  determined  for  attached 

fiow  on  the  flap  using  equation  (18).  Charts  showing  this  pressure  ratio  are  pre¬ 
sented  in  figure  18.  in  terms  of  Ax/6*HL.  and  the  flap  expansion  angle .  The 
xtrapoiation  factor  FiM  p  is  formed  for  any  location  on  the  flap  by  determining 
condidonsU.re  COrrespondi^  to  ^HL.  evaluated  at  wind-tunnel  or  fligM 

“e  hf-‘r  f"  *•  «*•  hinge  line  end  “XTT'’ 

the  present  analysis,  it  was  found  that  hfi*.  as  developed  in  reference  2  for  lir 

ipisiiissr 

tit  •sn *  Tk  ‘ir™  -  19 

d  "  i  T  °,  Figure  19  corresponds  to  free -stream  flow  temperatures  of  100 

pressibiUtv  \  *  iS  Sh°W"  in  figUre  20  in  ^ms  of  the  three  com¬ 

pressibility  factor-temperature  products  /7T\  nn  „  j 

script  w  denotes  evaluation  at  the  wall  The  nr  ’  a  )w.  and  (ZT)o-  The  sub- 

<*»  he-  in  .ho.™  .he  steShe.  b"ow.  P™“'“r<'  “»  ’™- 
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The  extrapolation  factor  FjMh  for  heat  transfer  is  defined  as 


cm/»T)g 


500 


T*  -  TOO  *"ir 


<h/hHL>FLT,  (Ax/6*)flt 
FiM,h  <h/hHL>WT,  (Ax/6*)wt 


The  ratio  (hAjjL)  ln  equation  (21)  is  evaluated  using  curves  presented  in 
this  report  in  exactly  the  same  way  as  the  pressure  ratio  of  equation  (20).  Charts 
showing  the  heat  transfer  coefficient  ratio ,  h/hHL,  are  shown  in  figure  21.  In 
this  figure,  h/hjjL  Presented  in  terms  of  the  flap  expansion  angle,  at  a  10 
and  40  degrees  and  =  10  and  20. 

4.  INTERFERENCE  GEOMETRIES 

Maneuverable  reentry  vehicles  and  hypersonic  cruise  vehicles,  which  of 
necessity  are  designed  with  protruding  stabilizers  and  control  surfaces,  will  foster 
mutually  Interfering  flow  fields.  As  a  result,  some  portion  of  either  type  of  ve¬ 
hicle  may  experience  shock-wave  impingement,  or  separating  and  reattaching 
flows.  The  increased  pressure  associated  with  interfering  flow  fields  has  the 
tendency  to  reduce  the  boundary-layer  thickness  in  the  region  of  interference  and, 
consequently,  increase  the  local  convective  heat-transfer  rate.  Data  reported  in 
references  20  through  28  indicate  that  local  heating  rates  can  be  larger  by  an  order 
of  magnitude  or  more  compared  to  those  without  interference.  In  general,  the 
size  of  Interference  regions,  as  indicated  by  the  data,  is  small;  hence,  the  total 
heat  load  is  not  substantially  increased. 

An  interference  geometry  is  defined  as  a  body  which  generates  a  secondary 
shock  system  that  interacts  with  the  basic  shock  s>stem.  For  example,  a  swept 
cylinder-flat  plate  combination,  as  illustrated  below; 
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The  flat  plate  deflected  at  a  slight  angle  generates  a  secondary  shock  wave 
which  interacts  with  the  basic  shock  wave  generated  by  the  swept  cylinder.  The 
flat  plate  represents  the  slender  body  of  a  vehicle.  The  cylinder  represents  the 
leading  edge  of  a  wing,  a  stabilizer,  or  an  Inlet,  and  is  mounted  to  the  flat  plate 
at  a  sweep  angle  A.  The  cylinder  sweep  angle  is  always  measured  with  respect 
to  an  axis  perpendicular  to  the  free -stream  flow  direction. 


a.  Swept-Cylinder  Stagnation-Line  Interference  Heating 


The  data  extrapolation  method  developed  for  the  interference  geometry  is 
applicable  at  the  stagnation  line  of  the  swept  cylinder.  Stagnation -line  heat- 
transfer  data  that  has  been  measured  on  a  cylinder  in  a  disturbed  flow  field,  and 
normalized  using  swept -infinite -cylinder  theory,  is  extrapolated  from  wind-tunnel 
conditions  to  flight  conditions  using  the  following  format: 


NOTE:  Above  procedure  is  applicable 

to  both  laminar  and  turbulent  flows 


Extrapolated 

data 


FLT 
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The  wind-tunnel  data  is  expressed  in  the  form  (h/h*),  where  h  is  heat- 
transfer  coefficient  data  on  the  leading  edge  in  the  interference  region,  and  h  a 
is  the  stagnation-line  heat-transfer  coefficient  for  an  infinite  cylinder  of  the  same 
sweep  angle  and  radius.  The  factor  F;  is  used  to  extrapolate  the  data  from  wind- 
tunnel  to  flight  enthalpy  conditions  at  the  wind-tunnel  Mach  number.  The  factor 
FM  is  used  to  extrapolate  flight  values  from  the  wind-tunnel  Mach  number  to  the 
desired  flight  Mach  number.  The  procedure  is  identical  for  laminar  and  turbulent 
flows.  The  extrapolation  factor  Fi  is  defined  as: 


Fi 


^max^A^FLT 

_(hmax'/hAKvT_ 


M, 


WT 


(22) 


In  equation  (22),  the  maximum  heat-transfer  coefficient  in  the  vicinity  of  shock- 
wave  impingement,  hj^ax,  is  normalized  using  the  undisturbed  stagnation-line 
heat-transfer  coefficient  for  an  infinite  cylinder,  hA.  The  evaluation  procedure 
used  for  this  ratio  is  dependent  upon  the  flow-field  geometry.  The  first  method 
used  corresponds  to  a  flow  field  where  the  secondary  shock  wave  fostered  by  the 
cylinder  is  attached  at  the  intersection  of  the  cylinder  and  the  boundary -layer  re- 
attachment  location.  The  second  method  corresponds  to  a  flow  field  where  this 
secondary  shock  wave  is  detached.  The  flow  model  illustrated  in  the  following 
sketch  for  an  attached  shock  wave  was  first  postulated  in  reference  21. 
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For  this  flow  field,  the  ratio  (h^x/hy^)  is  defined  by  the  relationship 
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where  hmax  is  the  maximum  heat-transfer  coefficient  in  the  vicinity  of  shock  wave 
impingement,  hA  is  the  undisturbed  stagnation-line  heat-transfer  coefficient  for 
an  infinite  cylinder,  Pmax  is  the  stagnation-line  pressure  determined  by  the 
interference  method,  PA  is  the  stagnation-line  pressure  on  a  yawed  infinite  cylin¬ 
der,  and  the  subscript  L,  denotes  laminar  boundary-layer  flow.  The  subscript 
M„,  WT  in  equation  (22)  denotes  that  the  ratios  (hmax/hA)  for  both  flight  and 
wind  tunnel,  are  evaluated  for  a  free-stream  Mach  number  equal  to  that  in  the 
wind  tunnel.  Cylinder  stagnation-line  pressures  for  flow  with  interference  are 
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calculated  using  oblique  and  normal  shock  methods_as  required  for  the  flow  field 
illustrated.  The  term  [tsl>  (A  +  0)/TSLi  A]0'25  essentially  corrects  for  in¬ 
creases  in  the  effective  sweep  of  the  cylinder.  The  last  term  is  a  pressure 
gradient  term  tabulated  in  reference  29  in  terms  of  the  velocity  gradient  parameter, 
9.  In  this  analysis,  9  is  oqual  to  1.0  and  0.5  for  stagnation-line  and  stagnation- 
point  flows,  respectively.  Equation  (23a)  is  developed  in  reference  21. 


Turbulent  heat-transfer  coefficient  ratios  for  cylinder -plate  combinations 
with  attached  shock  waves  are  determined  using  the  relationship 


^hmaxj  ^  |pma xj'8  jgin  (A  +  a)j°'G  (23b) 

where  the  term  [sin  (A  +  a)/sin  a]0,6  corrects  for  the  velocity  component  normal 
to  the  cylinder,  and  the  subscript  T  denotes  turbulent  boundary-layer  flow.  Equa¬ 
tion  (23b)  is  presented  in  reference  21  and  is  a  modification  of  the  relationship 
presented  in  reference  30.  An  axial  pressure -gradient  term  which  appears  to  the 
one-fifth  power  in  the  latter  analysis  has  been  neglected. 

For  those  flow  conditions  where  the  cylinder  shock  wave  detaches  at  the  root 
of  the  cylinder,  the  flow-field  is  illustrated  by  the  following  sketch: 


For  this  situation,  the  laminar  and  turbulent  heat-transfer  coefficient  ratios 
are  defined  by  the  following  expressions: 


/hmax\  =  ^2  1  g'w,  Sp\ 

\  hA  jh  A1  (cos2  A)  \  g'w,  SLy 

/Nnax\  _  ^2  1 

\  hA  /  T  A1  (cos2  A) 


(24a) 


(24b) 
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Stag  ftownelTand  Z  tatae  ^ndoffdfstf  in  an  int°r- 

™8 — -  — -  52urssi! r-srit. 


where 


f  CkSL>  (r) 
<*SL> 


1  -  f 


<SL  =  Pm/p  SL 


and 


f  (Ksl)  =  ln  lg  (Ksl>  +  [k2  (ksl)  +  i]  1/2j 


1  -  kbi 

8  (KSL)  =  - 5L 


[kSL  (2  ~  KSL^J 


1/2 


(25a) 

(25b) 

(25c) 

(25d) 


«•«  srx-ifc”«r;  d°T*n  Tt  i  ‘,h“k  *»- 

cyu^r  station  Une,  »d  r  „  „„  dL“ iTJZS ‘he 

-r  • .r  :re  ‘r^ now 

expressed  1„  vt‘«„y  “i wuSfS “mPer,,,,re  r°ll°  “»  >» 

function  co82  A  for  hypersonic  free -stream  flow  mat®ly  reduces  to  016  sweeP  angle 

sSsSBS 

tunnels  is  shown  in  figure  24  The  dash  h  \i  or  “rbulent  n°w  ln  air-nitrogen 
detached  and  attactedSLk waJes Z? ^  ^  flows  with 

spectively.  Curves  arett  I  f  both  wind-tunnel  and  flight  conditions,  re¬ 
waves  in  wind  tunnels  and  attached  shoSk wa^^mS*' fetached  shock 
cussed  in  the  following  subsection.  The  cSves  in  fibres  J2  23  td?^ 
spend  to  a  total  enthalpy  ratio  I  =  0. 25,  where  ’  '  COrre- 


I  =  °-  WT 

io,  FLT 


(26) 
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The  extrapolation  factor  FM  is  defined  as 


(hr 


M 


(hn 


^A^M . 


FLT 


^A^M , 


WT 


(27) 


where  the  ratios  <hmaxAA)  are  obtained  from  figures  25  and  26.  The  subscripts 

fli’,h!  a  ’  !VT  lnd,1Cat°  that  the  heat  tra«sfer  ratios  are  evaluated  at  the 

flight  and  the  wind-tunnel  Mach  numbers,  respectively. 

tr  .n  Ir°  dcte™ine  absolute  heating  rates  in  flight,  the  swept-infinite -cylinder  heat- 
transfer  coefficients  must  be  obtained  from  a  separate  sou  rce.  Calculations  for 
swept-infinite  cylinders  in  laminar  and  turbulent  flows  are  presented  in  reference  2. 

b.  Dissimilar  Flow  Fields 

„  n  Then,u  6  Sil°Ck  WaVe  fostered  by  ‘he  cylinder  detaches  from  the  root  of  the 
cylinder  ,  the  interference  flow  field  undergoes  a  change.  Because  the  theories 
for  predicting  maximum  heat  transfer  behind  attached  and  detached  shock  waves  do 
not  converge  upon  each  other  at  shock-wave  detachment,  a  discontinuity  occurs 
LefktTt  .  f, at  th°  Cyllnder  is  dependent  upon  the  cylinder  sweep  angle 
lip!  fanf  ,  '  thC  SUrface  Mach  number,  and  the  total  enthalpy  ’ 

Dependence  of  shock-wave  detachment  upon  total  enthalpy  fosters  a  situatio/which 

z  11,6  w“'ta"ei  ““  *■  «**»*> 

reached  When  the  cylinder  sweep  angle  is  small  enough  to  make  shtk  wa^  de 
2 7° show's  P°sslblc'.detac|'menc  may  occur  in  wind  tunnels  but  not  in  flight.  Figure 

Where Th  u  .  .P°lnts  whlch  ^present  the  approximate  boundaries  befween  regions 

where  this  situation  may  occur  for  flat-plate  angles  of  attack  between  10  and  30^ 
egrees  The  dashed  and  solid  lines  indicate  wind-tunnel  and  flight  conditions 
respectively.  Wind  tunnels  represented  by  this  chart  are  assumed  to  be  clpabie 
f  reaching  approximately  one-fourth  the  total  enthalpy  of  flight  at  the  same  Mach 

Z'.Jf  “»  °f  “•  ■*.«.  win/tuZrr.  M.ch" 

If  the  plate  (shock  generator)  is  at  a  10-degree  angle  of  attack  the  dividing  Doint 
between  attached  shock  waves  and  detached  shock  waves  in  the  wind-tunneUs  at  a 
cylinder  sweep  angle  of  approximately  31  degrees.  For  sweep  angles  g^ter  than 
degrees  the  shock  wave  is  attached,  and  at  sweep  angles  less  than  3gl  degrees 
the  shock  wave  is  detached.  In  flight  at  Mach  15,  the  dividing  poirHetweer 
attached  and  detached  shock  waves  occurs  at  a  sweep  angle  ejial  to  approximately 

•re  rZi'J ln  n'ght  "d 

wind-tunnel  and  flight  environments. 
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It  should  be  notec  that  the  boundary  for  the  geometry  in  helium  ground  facil¬ 
ities  at  a  -  30°  lies  between  the  boundaries  at  a  =  10°  and  20°.  This  occurs 
with  changing  a  due  to  the  combined  effect  of  changes  in  the  cylinder  deflection 
angle  (90°  -  A  -  Of)  and  the  flat -plate  surface  Mach  number,  Ms. 

5.  ROUGHNESS  ELEMENTS 

It  is  improbable  that  the  surface  of  a  reentry  vehicle,  or  hypersonic  cruise 
vehicle,  will  have  the  ideally  smooth  surfaces  typical  of  wind-tunnel  test  models. 
Manufacturing  tolerances,  load  deformation,  and  thermal  expansion  requirements 
will  cause  flow -field  disturbances.  Extrapolation  methods  have  been  developed 
for  two  types  of  surface  roughness  elements,  (1)  shallow  convex  waves,  and  (2)  butt 
joints  having  a  rectangular  cavity  below  the  contour  line  at  the  joint.  These  ele¬ 
ments  are  illustrated  by  the  sketch  below: 


Rectangular  Cavity 
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I  he  procedure  that  has  been  developed  for  extrapolating  shallow  convex  wave 
and  rectangular  groove  maximum  laminar  heating -rate  wind-tunnel  data  to  flight 
conditions  is  illustrated  by  the  following  sketch 


A  description  of  the  extrapolation  factor  F1M  and  methods  for  its  development 
follows. 


a.  Shallow  Convex  Waves 


In  order  to  extrapolate  the  maximum  laminar  heat -transfer  coefficients 
obtained  in  the  wind  tunnel  to  flight  conditions,  the  measured  values  must  first 
be  normalized  by  smooth  surface  heating  rates  at  the  same  distance.  The  ratio 
(h/ham)data  is  extrapolated  to  flight  by  the  extrapolation  factor  FiM  which  is 
defined  as 


F 


iM 


(A^/eism1  FLT.  R/fi* 

^^smW.R/e* 


(28) 


The  ratio  (Aq/qsm),  where  Aq  =  qmax  -  qsm,  in  equation  (28)  is  deter¬ 
mined  using  the  following  expression 


R_ 
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sm 


_ 1 _ 

2.5[.78  7  -  84  (iw/l0)] 
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9 

(W/R)  -y 

2  -  1~| 
e,  n  | 

•ny 

[-78  .  .84(iw/io,]| 

M 

2  * 

(29) 


where 


q  -  q 
-  _  max  Hsm 

q  “  6  <3°) 
sm 

The  derivation  of  equation  (29)  is  presented  in  reference  19. 

Curves  showing  solutions  for  q  on  a  sharp  flat  plate  in  air  are  presented  in 
figure  28  as  a  function  R/6*sm  and  [w/R  (Me  n2  -  l)°-5/M(1  n2].  These 
curves  are  presented  for  values  of  the  effective  specific  heat  ratio,  V  equal 
to  1.4  and  1.15.  The  parameter  R/68m*  is  evaluated  using  the  exprlssion 

R_  _  hR 

6*  (h6  *) 

sm  'sm 

where  h^^  Is  a  flat-plate  heat-transfer  coefficient  evaluated  at  a  distance  equal  to 
that  between  the  leading  edge  and  the  center  of  the  wave,  R  is  the  maximum 
height  of  the  wave,  and  the  product  hfl*  is  obtained  using  the  procedure  outlined 
in  Section  m-3.  Edge  Mach  number  in  the  parameter  [ W/R (M  2  -  l)°-5/M  2  ] 

is  evaluated  in  the  plane  normal  to  the  wave  using  the  expression"  e,n 

M  =  M  cos  (A)  no\ 

e,n  e  '  l/wave  w2> 

The  extrapolation  factor  FiM  for  waves  on  a  sharp  flat  plate  is  evaluated  using 
values  obtained  from  these  curves.  The  method  of  extrapolation  presented  here 
for  waves  takes  into  account  the  physical  characteristics  of  the  wave  with 
respect  to  the  predicted  depth  of  the  displacement  thickness  for  both  hot-wall 
and  cold-wall  conditions. 

Absolute  values  of  the  maximum  heating  rate  predicted  for  flight  conditions 
can  be  obtained  by  determining  hgm  using  reference  2  for  the  flight  condition 
of  interest. 

No  reliable  analytic  or  empirical  method  exists  for  extrapolating  turbulent 
wind-tunnel  data  to  flight  conditions  for  waves. 

b.  Rectangular  Grooves 


Extrapolation  of  maximum  laminar  heating-rate  measurements  obtained  in 
an  air-nitrogen  wind  tunnel  for  rectangular  grooves  is  accomplished  by  first 
normalizing  this  heating  rate  using  smooth-body  heat-transfer  coefficients 
evaluated  at  a  distance  equal  to  that  between  the  leading  edge  and  the  center  of 
the  groove.  The  extrapolation  factor  FiM  is  defined  as 


F 


iM 


^max^sm^FLT 

^maAmW 


(33) 
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the  rat'°  Nax^sm*  is  evaluated  using  the  empirical  correlation  expres¬ 
sion  for  a  rectangular  groove  that  was  first  presented  in  reference  18.  This 
correlation  is  expressed  as 


max 


=  1  +  (  0.035  + 


h  WPr 
am  w 


(34) 


equafon  (34).  W  is  the  width  of  the  cavity,  y  is  the  depth  of  the  cavity.  hBm 

COefflclent  on  a  smooth  body  evaluated  at  a  location  corre- 
pon  ing  to  the  mid-point  of  the  groove,  Prw  is  the  Prandtl  number  based  on 
wall  temperature  and  pressure,  Mw  is  the  viscosity  and  cp  w  is  the  specific 
heat  at  constant  pressure  evaluated  at  the  wall  temperature’. 

mi  fTeS  PureSen,ted  ln  flgure  2b  which  represent  solutions  for  equation 
(34)  in  terms  of  the  width  to  depth  ratio  of  the  groove,  W/y,  and  the  parameter 

h  WPr 
am  w 


Li  C_ 

W  P.w 


The  extrapo.ation  factor  F 1M  is  formed  by  evaluating  the  ratio  <h  /h  i 
for  values  corresponding  to  h  WPr  /u  r  nt  .  ,  'bmax'"sin) 

hw  ag  io  "smw^rw/4wcP,w  at  wind-tunnel  or  night  conditions. 

No  reliable  analytic  or  empirical  method  exists  for  extrapolating  turbulent 
wind-tunnel  data  to  night  conditions  for  grooves. 

6.  CHEMICAL  NONEQUILIBRIUM  EFFECTS  ON  EXTRAPOLATION  FACTORS 

_  f°r  the  free-stream  conditions  of  present  day  wind  tunnels  and  in  the  orbital 

the  eaT  Tm  k’nfC  energy  is  the  greater  part  of  the  total  energy  of 

the  gas.  When  this  gas  is  brought  to  rest  by  compression,  the  Row  energy  is 

converted  into  thermal  energy  resulting  in  a  considerable  increase  in  temperature 
molecular  excitation,  and  dissociation.  This  compression  may  be  a  result  of 
crossing  a  strong  shock  wave  generated  by  a  blunt  body,  or  of  viscous  inter¬ 
action  in  the  boundary  layer  of  a  slender  body.  These  two  types  of  compression 
are  illustrated  in  the  sketch  below: 
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Region  of  high  temperature 
effect* 


of  high 
temperature 
effect* 


:zrzz\zzizTon  due  * viscous  «»*»»««  -  * 

compared  to  their  nose  or  leading  edg^ Gamete!-08  ’  ^  long 

si iss=sE~: 

flow  field  are  difficult  to  assess  rltrorm.*^  J  “  °"  otherwise  well-defined 
less  confidence  In  order  tn  B,y  and  are  treated  analytically  with 

flow  fields  enveloping  the  geometrTeuLrtT^8  present  analysis.  the 

Sion  have  been  simplified  The  simnHrt  m  8  68Cr  in  the  prevlous  disc  us - 
treatment  of  therSem^  a"  -alytic 

developed  for  these  elements  How  e.Cts  on  tlle  extrapolation  factors 

polation  factors  mly  te  affect'ed  bv  r"'  d‘SCUSsion  of  ho-  the  ext  ra¬ 
in  this  section.  X  ther™che"^al  nonequilibrium  is  included 


a*  Blunt  Bodies 

surfaces  deflectedhLto  the  Stream*  ™*  *"d  g~,  -ntrol 

sharp  delta  wings  at  angles  of  attack  gre^teJ thfnatoS  Sd  C°mbina.tion8’  and 
for  air-nitrogen  wind-tunnel  conditions  than  it  is at  Htelrt  coIXf68’  ,8/reater 
result,  heating  rate  predictions  for  fiic-ht  k  a  ^  ‘  ght  condltions.  As  a 
sented  in  this  report  mav  he  hitrv,  nf  ’  niSe  on  extrapolation  iactors  pre¬ 
conditions.  V  ^  high6r  tha"  WiI1  actual,y  be  experienced  at  night 


co  abov'’  ,omr  ^ 

occur.  Near  the  .uZ^  to 

0r  ^  ™dl“'  "  •»=  «».  «™/be 
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than  a  nonequilibrium  state.  Actually,  in  this  region,  the  gas  is  rapidly 

— -fio  eqU‘libriUni  state  1,6031186  of  increase  in  collision  frequency. 

8  the  n°w  pr°gresses  away  from  the  stagnation  region,  it  undergoes  a 
apid  expansion  and  could  pass  through  a  series  of  quasi-frozen  states  In  a 

SsTdTrri’biee  sr  r rate  is  so  smau  z  the 

•i  u,  negligible.  It  is  shown  in  reference  31  that  the  reduction  in  the 
availaWe  energy  can  result  in  decreases  on  the  order  of  60%  in  heat  transfer  to 

y  compared  to  that  which  would  occur  for  equilibrium  now-field  conditions. 

b.  Slender  Bodies 

In  air-nitrogen  wind  tunnels  (excluding  combustion-driven  shock  tunnels) 

sharp  delta  wing  wind-tunnel  data  are  negligible. 

Sharp  delta  wings  at  low  angles  of  attack  and  at  night  conditions  will  cause 
ow  compression  through  viscous  interaction  in  the  boundary  layer,  since  the 
majority  of  the  kinefc  now  energy  of  the  free-stream  has  Jen  retaineTac«.ss 
the  relatively  weak  shock  wave  generated  by  this  geometry  As  the  now 

F°r  sharp  cones,  chemical  nonequilibrium  nonsimilar  boundary-layer  oro- 
gram  calculations  indicate  that  temperature  gradients  at  the  wall  are  less  than 
would  be  expected  for  gases  in  chemical  equilibrium  as  the  now  progresses  from 
e  apex.  This  reduction  in  temperature  gradient  at  the  wall  fosters  changes  in 
the  specie  concentration  profile,  which  in  turn  leads  to  more  specie  difftS  4 

reference  IT  Tthat  rfjj"  °"  8l6nder  ***  ^  transfer  is  bussed  in  ' 

1,17  ,  *  h  eport’  gas  specie  dlffusi°n  is  shown  to  be  of  relatively 

greater  importance  in  the  slender  body  case  than  in  the  blunt  body  case  This 

Decreases  up* to  ^  ^  U  mUch  thlnner  near  0,6  wall. 

mtehtX  *  Tu  beat  transfer  due  to  gas  specie  diffusion  on  slender  bodies 
might  be  expected  based  on  the  results  of  reference  35.  This  result  does  not 
agree  with  the  results  shown  in  reference  32  where  chemical  nonequihbrium 

fn  tha,  "  gas  8Pect6  dlffuslon  was  shown  to  increase  heating  slightly.  However 
n  hat  report,  the  slight  increase  shown  was  offset  by  the  decrease  resulting  ’ 
rom  the  decrease  in  the  temperature  gradient  at  the  wall  with  increasing  stream 


3  See  Appendix  B. 

4  A  more  detailed  discussion  is  contained  in  references  32,  33  aw*  34 
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c.  Effects  of  Catalytic  Wall 


By  definition,  a  catalytic  wall  fosters  recombination  of  dissociated  gas 
molecules  in  proximity  to  the  wall,  thus  releasing  the  energy  absorbed  in  disso¬ 
ciation.  Heat  tra.isfer  to  surfaces  enveloped  in  a  flow  field  that  is  dissociated 
is  therefore  a  function  of  wall  catalycity. 

For  the  slender  bodies  considered  in  this  report,  the  effect  of  wall  cataly¬ 
city  on  heat  transfer  is  always  small,  since  the  degree  of  dissociation  is  small 
for  both  flight  and  wind-tunnel  conditions.  On  the  other  hand,  a  significant 
increase  in  nonequilibrium  heat  transfer  to  blunt  bodies  can  occur  because  flow 
compressed  by  a  strong  shock  wave  can  cause  a  high  degree  of  molecular  disso¬ 
ciation.  These  molecules  tend  to  recombine  in  proximity  to  a  catalytic  wall. 

It  is  indicated  in  reference  32  that  for  blunt  bodies  in  a  highly  dissociated  flow 
field,  a  noncatalytic  wall  can  reduce  heat  transfer  at  the  stagnation  point  by 
as  much  as  75%  in  comparison  to  a  catalytic  wall. 

Increases  in  heating  that  can  be  attributed  to  catalycity  has  been  accounted 
for  in  the  extrapolation  factors  presented  in  this  report.  All  calculations  were 
made  for  a  fully  catalytic  wall,  and  apply  to  wind-tunnel  data  measured  on 
models  that  have  highly  catalytic  surfaces. 
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SECTION  rv 


CONCLUDING  REMARKS 


One  purpose  of  this  analysis  was  to  evaluate  the  capability  of  operational 
ground  facilities  to  simultaneously  simulate  the  free-stream  Mach  number 
Reynolds  number,  Knudsen  number,  Lewis  number,  and  enthalpy  (velocity) 
that  will  occur  during  hypersonic  reentry.  It  was  found  that  only  free-stream 
Mach  number,  free-stream  Reynolds  number,  free-stream  Knudsen  number, 
and  boundary-layer  edge  Lewis  number  are  simulated  adequately  in  present- 
day  facilities.  Flight  condition  boundary-layer  edge  Mach  number  and  Reynolds 
number,  two  parameters  that  must  be  simulated  exactly  to  insure  similar 
boundary-layer  growth  at  ground-test  conditions,  are  not  simulated  well  in 
present-day  ground  facilities.  The  largest  deficiency,  however,  is  due  to  the 
failure  to  adequately  simulate  high  Mach  number  flight  enthalpy  conditions. 
Present-day  facilities  large  enough  for  development  testing  have  the  capability 
to  run  at  only  about  half  the  velocity  occurring  during  the  hottest  portion  of 
reentry. 


The  primary  purpose  of  this  analysis  was  to  develop  extrapolation  factors 
using  existing  analytic  methods  that  account  for  any  simulation  deficiencies 
occurring  in  present-day  ground  facilities.  These  factors  were  to  be  used  in 
the  application  of  ground-test  data  to  reentry  vehicle  design  and  were  to  be 
developed  for  four  geometric  elements  that  are  representative  of  portions  of 
a  reentry  vehicle.  The  elements  are  a  sharp  unyawed  delta  wing,  deflected 
control  surfaces  In  compression  and  expansion,  shallow  convex  waves  and  rec¬ 
tangular  grooves,  and  a  yawed  circular  cylinder  intersecting  a  sharp  flat  plate. 

Extrapolation  factors  have  been  developed  for  these  four  elements  that  are 
based  on  existing  analytic  methods.  However,  It  was  found  during  the  present 
analysis  that  the  best  analytic  methods  for  these  elements  are  only  approximate 
and  that  for  the  groove  the  best  method  is  empirical.  In  addition,  existing 
methods  generally  neglected  real -gas  effects. 


x  .u  newiflow-fleld  method  was  developed  for  a  thin-sharp-unyawed  delta  wing, 
hi  this  method,  flow-field  gas  properties  in  proximity  to  the  centerline  of  the 
delta  wing  are  computed  using  a  new  similarity  parameter  that  accounts  for 
both  real-gas  and  leading-edge  sweep-angle  effects.  Away  from  the  centerline 
at  angles  of  attack  greater  than  about  50  degrees,  the  present  delta-wing  method 
is  dependent  upon  correlations  of  the  approximate  analytic  method  presented  in 
reference  36. 


The  delta-wing  flow-field  method  presented  here  postulates  that  a  parting 
line  exists  on  a  sharp  delta  wing  at  angles  of  attack  between  about  30  to  60 
degrees  that  accounts  for  both  in-flow  and  out-flow  occurring  on  the  surface 
simultaneously.  The  present  analysis  Indicates  that  maximum  span  wise  heating 
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°£  the  difference  increasing  with  increasing  deflection  angles.  The  distributions 
also  Indicate  that  differences  between  surface  pressures  on  the  expansion  surface 
obtained  in  a  wind  tunnel  and  those  obtained  in  flight  may  be  greater  than  100% 
over  some  portion  of  the  flap.  There  are  two  reasons  for  this:  first,  the  total 
pressure  drop  for  a  given  deflection  angle  is  greater  in  flight  and  second,  the 
boundary  layer  at  the  hinge  line  in  flight  is  much  thicker  than  it  is  in  a  wind 
tunnel.  Additional  analytic  analysis  is  required  to  more  precisely  determine 
the  effects  of  the  boundary  layer  on  the  pressure  distribution  and  conversely 
the  effects  of  the  pressure  distribution  on  the  boundary-layer  development. 

Further  development  of  the  analytic  method  for  predicting  heat-transfer  to 
Ulff  CyIindG1S  Intersecting  a  sharp  flat  plate  is  required.  The  method 
used  in  this  document  predicts  a  large  discontinuity  in  pressure  and  heat-transfer 
when  the  shock  wave  fostered  by  the  cylinder  detaches  from  the  loci  of  intersec¬ 
tion  with  the  plate.  This  is  reasonable,  but  methods  are  required  to  establish 
how  to  interpolate  between  pressures  and  heating  rates  that  occur  when  the  shock 
wave  is  attached  at  wind-tunnel  conditions  and  detached  at  flight  conditions. 

Differences  between  wind-tunnel  and  flight  heating  rates  for  the  cylinder- 
plate  combination,  for  surface  waves  and  grooves,  and  for  deflected  flaps  have 
been  found  to  be  extremely  geometry  dependent.  Extrapolation  factors  presented 
in  this  report  vary  from  unity  to  values  exceeding  50%  for  some  conditions. 

In  general,  the  methods  chosen  for  development  of  the  extrapolation  factors 
for  these  four  geometric  elements  reflect  the  present  state-of-the-art  As 
more  precise,  or  better  substantiated,  analytic  methods  become  available  new 
extrapolation  factors  should  be  developed.  The  procedures  developed  in  the 
present  analysis  and  presented  in  this  report  should  be  applicable  at  that  time. 
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Figure  h  VELOCITY  SIMULATION  IN  GROUND  FACILITIES 


Figure  2:  ENERGY  SIMULATION  IN  GROUND  FACILITIES 
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Figure  3:  FR£E-$TR£AM  REYNOLDS  NUMBER  SIMULATION  IN  GROUND  FACIHJfES 
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Figure  4:  REYNOLDS  NUMBER  AND  INVISCID  MACH  NUMBER 
COMPARISON  ON  A  SHARP  FLAT  PLATE 
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Figure  5:  REAL-GAS  EFFECTS  ON  REYNOLDS  NUMBER  FOR  WEDGES  AND 

SHARP  UNYAWED  CONES 
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Figure  5?  Concluded 


Figure  6:  ORBITAL  REENTRY  CORRIDOR  FOR  LIFTING  VEHICLES 
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Figure  h  SHARP  delta  WING  LAMINAR  EXTRAPOLATION  FACTORS  FDR 


Figure  8:  SHARP  DELTA  WING  LAMINAR  EXTRAPOLATION 
FACTORS  FOR  HELIUM  GROUND  FACILITIES 
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Figure 
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Figure  10:  SHARP  DELTA  WING  CENTERLINE  REYNOLDS  NUMBER  DISTRIBUTION 


49 


Continued 


c)A=  80° 

Figure  10:  Concluded 


Figure  Ih  REFERENCE  REYNOLDS  NUMBER  FOR  AIR  AND 
NITROGEN  GROUND  FACILITIES 
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Figure  12:  Concluded 
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Figure  13:  REYNOLDS  NUMBER  FUNCTION 


On  25 
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a>  A  -  70°  b)A=  80o 

Figure  14:  SHARP  DELTA  WING  SPANWISE  LAMINAR  EXTRAPOLATION  FACTORS 
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Figure  15:  SHARP  DELTA  WING  CENTERLINE  SONIC  FLOW  CONDITIONS 
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Figure  16:  MAXIMUM  HEATING  RATE  INCREASES  FOR  POSITIVELY  DEFLECTED  CONTROL  SURFACES 

IN  AIR  AND  NITROGEN  FLOWS 
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<2*/fV 

Figure  17:  MAXIMUM  HEATING  RATE  INCREASES  FOR  POSITIVELY  DEFLECTED 
CONTROL  SURFACES  IN  HELIUM  GROUND  FACILITIES 


T„  =  100  and  400  °R 


Figure  18:  PRESSURE  DISTRIBUTION  ON  NEGATIVELY  DEFLECTED  CONTROL 

SURFACES  IN  AIR  AND  NITROGEN  FLOWS 
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Figure  20:  DISPLACEMENT-THICKNESS  HEATING  PARAMETERS 
FOR  SHARP  FLAT  PLATES 


Figure  21 
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a  )a=  io°  b)a=  200  c)  a  =  s 

Figure  23:  INTERFERENCE  GEOMETRY  LAMINAR  EXTRAPOLATION  FACTORS  FOR 

HELIUM  GROUND  FACILITIES 
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a)  a=  io  b)  a=  20°  c)a  =  30° 

Figure  24:  INTERFERENCE  GEOMETRY  TURBULENT  EXTRAPOLATION  FACTORS  FOR  AIR  AND  NITROGEN 

GROUND  FACILITIES 


ill  EEL*.  ^ 


:::: 


^iJiEElDllima 


4ffimn 


:V;iyiiJljlu_Ji| 


Orbital  reentry 

Air  and  nitrogen  ground  facilities 


boundaries  of 
dissimilar  flow  fields 
a  =  10° 


30 

A  ,  degree* 

HJ 


0.25 


<'o)wr/(i0VLT  ~ 

- Helium 


Figure  27:  BOUNDARIES  BETWEEN  REGIONS  OF  ATTACHED  AND 

DETACHED  SHOCK  WAVES 
FOR  INTERFERENCE  GEOMETRIES 
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Figure  28:  SHALLOW  SURFACE  WAVE  LAMINAR  HEATING 
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c)  iw/i0  =  0.5 

Figure  28:  Continued 
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Figure  29:  MAXIMUM  LAMINAR  HEATING  NEAR  RECTANGULAR 

GROOVES 
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APPENDIX  A 


HEAT  TRANSFER  PREDICTION  METHOD 
FOR  SHARP  UNYAWED  DELTA  WINGS 


The  heat  transfer  analysis  for  the  sharp  unyawed  delta  wing  is  made  using 
the  Prti.  momentum  integral  method  described  in  references  2  and  18.  A  solu¬ 
tion  is  obtained  using  this  method  by  numerically  evaluating  integral  functions 
along  a  streamline.  This  solution  is  dependent  upon  the  thermodynamic  gas 
properties,  the  wall  temperature,  the  velocity,  and  the  cross-flow  pressure 
gradient  along  the  streamline.  When  any  one  of  these  quantities  varies  with 
the  streamline  distance,  the  solution  is  also  dependent  upon  the  upstream  his¬ 
tory  of  the  boundary  layer  flow. 

Along  streamlines  with  constant  flow-field  properties  and  with  a  finite,  but 
constant,  cross-flow  pressure  gradient,  the  integral  functions  are  evaluated 
explicitly,  and  the  entire  integral  solution  is  reduced  to  a  point  solution.  These 
requirements  are  satisfied  at  the  centerline  of  a  sharp  unyawed  delta  wing  at 
angle  of  attack. 


In  the  following  discussion,  the  explicit  form  of  the  solution  is  developed 
from  the  integral  form  for  the  centerline  of  a  sharp  delta  wing.  The  methods 
for  calculating  spanwise  heating  at  both  low  and  high  angles  of  attack  are 
described  in  the  last  section. 


1.  METHOD  AT  THE  CENTERLINE 


The  general  expression  for  evaluation  of  laminar  or  turbulent  heat  transfer 
to  an  arbitrary  geometry  is 


_ a. 

i 

aw 


i 

w 


=  H  = 


^  M0  Q 

0.645 

C 


(Al) 


where  <^/iaw  -  iw  is  the  enthalpy  heat -transfer  coefficient,  2  is  a  diffusion 
parameter  reflecting  the  effect  of  species  diffusion  on  heat  transfer  in  disso¬ 
ciated  How,  is  the  reference  stagnation  viscosity,  ar  is  the  reference  partial 
Prandtl  number,  QL  is  a  laminar  streamwise  pressure  gradient  profile  parame¬ 
ter,  Fx  and  xeq  L  are  equivalent  distance  parameters,  and  f{R_)  is  a  Reynolds 
number  function  expressed  ns 

fey  =  0.332  {rt  (A2a) 

for  laminar  flow,  and 
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(A2b) 


0.185  P, 


f(RjJ 


[log10  <Rr  +  3000)]  2> 
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for  turbulent  flow.  Rr  is  the  local  reference  Reynolds  number  defined  as 

0  M  u  x 

R  =  — — — - — - — £3j_L 

r  2  _  2 

H  F 


(A3) 


where Pj.tij-  is  the  reference  density- viscosity  product,  and  ug  is  the  boundary 
layer  edge  velocity. 


In  the  reduction  of  the  integral  solution  to  a  point  solution,  the  parameters 
defined  by  equations  (A2)  and  (A3)  are  not  affected.  The  profile  parameter,  QL, 
and  distance  parameters,  Fx  and  x  however,  are  integral  expressions 

which  are  modified  in  order  that  they’can  be  solved  explicitly.  In  the  integral 
solution,  the  profile  and  distance  parameters  are  expressed  as 


=JL 


3/10 


x  t  =  JT 
eq,  L  L 


-4/10 


F  =  J. 
x  L 


17/20 


[beq,  l] 
[beq,T/beq,L,] 


(A  4a) 
(A4b) 

(A4c) 


where  J  L  is  a  streamwise  pressure  gradient  parameter  of  relatively  complex 
form.  For  the  centerline  of  the  sharp  delta  wing,  JL  assumes  the  value  unity 
and  for  that  reason  its  definition  is  not  shown.  The  expression  for  JL  is  pre¬ 
sented  in  reference  18.  The  parameters  b„_  T  and  b  „  ~  are  expressed  as 

c4  *  ■Lj  GC|  y  1 


_  fX  /  E  \2  (A5a) 

beq,  L  =  [X]  /  PrUr  “e(rf  L)  d  x 


eq,  T 


"  M  I  »r»r“e(« 


M5/4 


d  x 


(A  5b) 


where  r  and  f  are  divergence  parameters  defined  later  in  this  Appendix,  EL  and 
E>p  are  the  laminar  and  turbulent  spanwise  pressure  gradient  parameters, 
respectively,  and  dx  is  d(x/xf).  The  subscript,  i,  refers  to  the  local  point  of 
interest.  The  bar  above  the  integrand  in  the  expressions  for  beg  L  and  be  „ 
Indicates  that  the  integrand  is  normalized  with  respect  to  the  value  at  x  =1tj! 
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At  the  centerline  of  a  sharp  delta  wing  the  streamwise  pressure  gradient  is 
zero,  and  the  thermodynamic  gas  properties  are  constant.  Hence, 

JT  s  i  (A 6) 

L 


(b  /x) .  = 
eq ,  L  1 


,  ■  / 


^eq.T^i  *  jf  (rfET)"  '  d  5 

It  remains  to  express  the  divergence  parameters  r  and  f  in  the  integrands 
in  terms  of  the  variable  x  and  to  evaluate  these  integrals. 

Evaluation  of  r  and  f  is  made  by  assuming  different  functions  of  x  that 
depend  on  the  shape  of  the  streamline  patterns.  In  the  vicinity  of  the  centerline, 
the  streamlines  take  the  form 

en 

T?  —  4 

The  exponent,  n,  is  a  crossflow  pressure  gradient  divergence  parameter  evalu¬ 
ated  at  the  centerline.  It  is  defined  in  Appendix  B.  Possible  streamline  pat¬ 
terns  are  illustrated  in  the  following  sketch: 


- 


When  n  >  1,  the  delta  wing  streamline  pattern  is  identical  to  the  streamline 
pattern  for  a  sharp  yawed  cone  and  the  parameters  f  and  r  become 


7  —  n-1 

f  =  x 


n  >  1 


For  values  of  n  <  1,  all  delta-wing  streamlines,  except  those  coinciding  with 
the  centerline,  originate  at  the  leading  edge  of  the  wing  rather  than  at  the  apex. 
These  streamlines,  therefore,  pass  through  regions  of  varying  pressure,  and 
strictly  speaking,  cannot  be  treated  using  methods  which  assume  that  a  zero 
streamwise  pressure  gradient  exists.  Sufficiently  far  downstream  and  near  the 
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centerline,  however,  this  effect  can  be  neglected.  All  necessary  flow-field 
requirements  are  satisfied  by  assigning  the  following  relationships  to  f  and  f . 


f  =  x  n(n-1) 


_  n(2-n) 
r  =  x 


n  <  1 


With  the  relationships  shown  in  (A8)  and  (A9),  equation  (A 7)  reduces  to 

l-l 


nil 
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^eq.T7^  =  (Xeq,T/x)  =  1  +  4  1  +  <n  -  D  ET 


(b  /x)  =  (x  /x) 
eq,  L  eq,  L  ' 


n  <  1 


=  |  1  +  2n^2  -  n  +  (n  -  1)  ELjJ 
^eq.T7^  =  <Xeq,T/x>  =  f1  +  \  n  [2  ’  n  +  (n  -  D  E 


(A  9) 


(A  10a) 

(A  10b) 

-1 

(A  10c) 


r 


(A10d) 


Ej^  and  E-j.  are  functions  of  thermodynamic  variables  only  and  can  be  evalu¬ 
ated  at  any  point  of  interest  by  the  relationships  shown  below. 
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(A  lib) 

(A12a) 

(A  12b) 

(A13a) 

(A13b) 

(A  14a) 

(A14b) 


(A15a) 
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For  n  >  1, 
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k  =  f  (n  -  1) 


and  for  n  <  1. 
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k  =  3  n(n  -  1) 


(A15b) 

(A  16a) 

(A  16b) 

(A16c) 

(A  17a) 
(A17b) 
(A  17c) 


Heat-transfer  equations  for  laminar  and  turbulent  flow  are  expressed  in 
terms  ot  the  above  parameters.  The  laminar  heat-transfer  expression  is  obtained 
by  combining  equations  (Al),  (A2a),  and  (A3).  With  a  slight  rearrangement  this 
reduces  to  the  form 


_  0.332  X 

hL  -  0.645 

r 


p  u  u 
r  r  e 


0.5 


(A18) 


where  x  is  the  centerline  distance  measured  from  the  apex  and  (xeq>  j^/x)  is 
obtained  using  equation  (A10). 


In  the  absence  of  cro3s-flow  pressure  gradients,  n=0  and  (Xgq^/x)  = 
1.  It  is  seen  that  the  expression  for  heat  transfer  at  the  centerline  of  a  sharp 
delta  wing  is  composed  of  an  expression  without  the  cross-flow  pressure  grad¬ 
ient  effect  which  is  multiplied  by  the  factor  (xeqi  l/x)”®’ J.  This  distance  fac¬ 
tor  is  the  only  factor  which  incorporates  the  cross-flow  pressure  gradient 
effect  and  it  is  evaluated  using  the  streamline  divergence  parameter,  n. 
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The  turbulent  heat -transfer  expression  is  based  upon  a  relationship  be¬ 
tween  the  laminar  and  turbulent  forms  of  the  Reynolds  number  function  expressed 
by  equation  (A2a)  and  (A2b).  Laminar  heat-transfer  coefficients  obtained  using 
equation  (A18)  are  modified  to  obtain  the  turbulent  heat-transfer  coefficient,  by 
the  following  expression 

0.185/ir 

”T  =  "L  0.332  flog10  (Rr+r3000)]2-584 


(A  19) 


where  Rr  is  evaluated  using  equation  (A3),  Equation  (A3)  after  slight  modifi¬ 
cation,  can  be  expressed  as 
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_r~r 

n  2 


F  2 
x 


(A20) 


where  F 

x 


is  evaluated  by  the  expression 


(Xeq,T/*) 

<xeq,L/x> 


1 

3 


(A21) 


The  reference  Reynolds  number  is  also  dependent  on  the  streamline  divergence 
parameter  n. 


2.  METHOD  AT  SPANWISE  LOCATION 


a.  At  Low  Angles  of  Attack 

The  spanwise  heating  distribution  at  low  angles  of  attack  can  be  expressed 
in  terms  of  centerline  heating  using  geometric  considerations  alone.  Thermo¬ 
dynamic  gas  properties  and  the  velocity  are  assumed  to  be  constant  on  the  lower 
surface  so  that  the  heat-transfer  decay  is  a  function  of  distance  only. 

Consider  the  streamline  "AB"  in  the  following  sketch: 


The  point  on  the  centerline  corresponding  to  the  point  ”B"  is  at  "D"  so  that  the 
comparable  centerline  distance  is  CD.  Because  all  properties  have  been  assumed 
constant,  the  heat  transfer  relationship  at  "B"  which  intersects  the  ray  angle  <f> 
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in  terms  of  that  at  "D"  which  lies  along  the  ray  angle  </>  0  can  be  written  as: 


h0=O 


(A22) 


where  m  has  the  value  0.5  in  laminar  flow  and  0.2  in  turbulent  flow  as  in  classi¬ 
cal  flat  plate  theory. 


The  distance  'AB  1  can  be  calculated  once  the  spanwise  distribution  of  the 
outflow  angle,  0  ,  is  known.  The  outflow  angle  distribution  is  assumed  to  vary 
according  to  the  relationship  0/8  *  =  (8//3!r .  At  low  angles  of  attack,  the  span- 
wise  heat  transfer  distributions  are  relatively  insensitive  to  the  value  of  the  ex¬ 
ponent  £  ,  in  this  expression  within  the  limits  0.5  i  C  <  2.  This  is  illustra¬ 
ted  in  the  sketch  below  for  three  values  of  C,  at  =  20  and  a  =  10  degrees 

and  20  degrees. 


0 
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The  differences  between  the  distributions  computed  for  the  range  of  C  between 
0.5  and  2  are  typical  for  the  range  of  sweep  angles  and  Mach  numbers  considered 
in  this  study.  In  extrapolation  factor  computations,  the  same  value  of  the  expo¬ 
nent  $  was  used  for  calculation  of  spanwise  heat-transfer  at  wind-tunnel  and 
flight  conditions 

b.  At  High  Angles  of  Attack 

The  PrPr  momentum  integral  method  is  applied  in  the  analysis  of  spanwise 
heat  transfer  at  high  angles  of  attack.  The  solution  is  obtained  numerically  by 
forward  integration  because  the  boundary-layer  edge  properties  and  the  cross- 
flow  pressure  gradient  are  no  longer  constant  along  the  streamline.  The  computa¬ 
tion  is  made  in  two  parts..  First,  a  streamline  is  calculated  on  the  lower  surface 
and  the  thermodynamic  gas  properties,  the  velocity,  and  the  cross-flow  parameters 


along  this  streamline  are  determined.  Second,  these  variables  are  used  to  compute 
the  heat-transfer  distribution  along  the  streamline  using  the  integral  form  of  the 
prHr  heat -transfer  method.  The  streamline  is  calculated  from  the  distribu¬ 
tions*^  of  the  local  outflow  angle  shown  In  Appendix  C.  In  most  cases,  a  number 
of  streamlines  must  be  calculated  in  order  to  obtain  solutions  for  all  spanwise 
locations. 
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APPENDIX  B 


METHOD  OF  FLOW  FIELD  SOLUTION  IN  PROXIMITY 
TO  THE  CENTERLINE  FOR  SHARP  UNYAWED  DELTA  WINGS 


Previous  methods  for  computing  thermodynamic  gas  properties  on  sharp 
delta  wings  employed  Newtonian  or  modified-Newtonian  techniques  which  neither 
differentiated  between  real  and  ideal -gas  flows  nor  accounted  for  variations  In 
leading-edge  sweep  angles. 

The  method  presented  in  this  appendix  was  developed  to  compute  thermo¬ 
dynamic  gas  properties  at  the  surface  in  proximity  to  the  centerline  of  a  sharp 
unyawed  delta  wing  at  angle  of  attack.  It  is  strictly  applicable  to  delta  wings 
having  leading-edge  sweep  angles  in  the  range  70  to  80  degrees  in  ideal  or  real- 
gas  flows  and  thermochemical  equilibrium.  This  method  makes  use  of  a  simi¬ 
larity  parameter  which  in  essence  permits  the  calculation  of  the  gas  properties 
near  the  centerline  without  using  the  method  of  characteristics.  The  value  of 
the  similarity  parameter  is  based  upon  correlations  of  analytic  and  experimental 
results. 


Additionally,  a  technique  is  developed  for  evaluating  cross-flow  pressure 
gradient  effects  at  the  centerline. 

1.  DEVELOPMENT  OF  THE  SIMILARITY  PARAMETER 


Flow-field  solutions  are  normally  difficult  to  obtain  for  complex  geometries 
without  some  idealization  or  mathematical  simplification.  Quite  often  a  com¬ 
plex  flow  field  can  be  separated  into  parts,  where  each  part  can  be  analyzed 
using  an  idealized  flow-field  model.  For  example,  the  unyawed  sharp-cone 
flow  field  can  be  separated  into  two  parts,  a  two-dimensional  shock  wave  and 
an  isentropic  shock  layer  between  the  shock  and  the  body  surface.  Normally, 
solutions  for  sharp-cone  flows  require  the  use  of  the  method  of  characteristics 
and  standard  shock  crossing  relationships  (see  for  example  references  40  and 
41).  Both  flow-field  parts  must  be  solved  simultaneously.  If,  however,  some 
information  about  each  of  the  parts  were  known  beforehand,  a  solution  could  be 
obtained  in  a  simpler  manner. 


Upon  examination  of  existing  solutions  across  shock  waves  generated  by 
wedges  and  sharp  unyawed  cones,  it  was  observed  that  the  parameter  denoted 
by  S  as  shown  below  assumes  a  unique  numerical  value  for  each  geometry. 


„  _  ^2_  tan  (8  -  op 
tan  0 


(Bl) 
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where  P ,  6,  and  a  are  the  density,  the  shock  wave  angle,  and  the  deflection 

angle, respectively.  These  parameters  are  illustrated  for  wedges  and  sharp 
cones  in  the  sketches  below: 


Wedge 

The  subscript  2  denotes  conditions  immediately  downstream  of  the  shock  wave. 
For  wedge  flows,  S  always  assumes  a  value  of  unity  because  the  expression 
for  it  is  a  form  on  the  continuity  equation  for  steady  two-dimensional  oblique 
shock  flows.  Values  of  S  calculated  from  exact  solutions  for  sharp  cones  are 
shown  in  figure  B1  as  a  function  of  the  Mach  number  normal  to  the  body  surface. 
In  the  regime  where  the  flow  normal  to  the  body  is  supersonic,  it  is  seen  that 
the  value  of  S  is  approximately  0.  5. 

The  similarity  parameter,  S,  when  once  defined  for  a  particular  geometric 
element  can  be  used  to  predict  some  surface  gas  property.  In  the  present  analy¬ 
sis  surface  pressure  is  correlated  to  S.  This  correlation  is  discussed  later 
in  this  section. 

The  fact  that  S  is  nearly  constant  and  independent  of  real-gas  effects  for 
both  wedges  and  sharp  cones  suggests  that  it  may  also  be  Independent  of  real- 
gas  effects  for  a  sharp  delta-wing  flow  field  in  proximity  to  the  centerline  where 
the  shock  wave  has  zero  curvature.  For  sharp  delta  wings  at  high  angles  of 
attack,  outflow  occurs  over  the  leading  edges  and  affects  the  flow  field  at  the 
centerline.  The  volume  of  outflow  is  dependent  upon  real -gas  effects  so  that  a 
change  in  out-flow  can  change  the  value  of  S  at  the  centerline.  It  will  be  shown 
that  while  the  real-gas  effect  on  S  is  significant  at  high  angles  of  attack,  the 
resultant  effect  on  thermodynamic  gas  properties  at  the  surface  is  not.  Despite 
this  relatively  small  resultant  effect,  a  real -gas  correction  to  S  has  been  used 
in  the  real -gas  analysis. 

*•  Similarity  Parameter  for  Sharp  Delta  Wings 


Variation  of  the  similarity  parameter  for  sharp  delta  wings  is  shown  in 
figure  B2.  At  high  angles  of  attack,  the  values  for  the  similarity  parameter 
were  correlated  to  flow-field  solutions  obtained  using  the  method  of  reference 
36.  At  low  angles  of  attack,  values  for  S  were  obtained  through  correlations 
of  surface  pressure  data  (references  42  and  43)  and  schlieren  photographs  of 
shock  shape  (reference  44).  Because  pressure  data  generally  is  accurate  to 
about  ±5%,  some  scatter  in  the  values  calculated  for  S  can  be  expected.  The 
values  calculated  for  S  are  very  sensitive  to  pressure.  This  is  shown  in 
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figure  B2  by  the  band  that  Indicates  the  variation  In  S  due  to  a  5 %  change  in  sur¬ 
face  pressure.  Nevertheless,  this  sensitivity  indicates  that  a  high  degree  of 
accuracy  in  S  is  not  necessary,  since  a  sizable  error  in  S  will  only  produce  a 
moderate  error  in  flow-field  and  surface  properties. 

At  low  angles  of  attack,  S  approaches  unity  (i.  e. ,  the  conical  shock  wave 
and  flow  field  begin  to  resemble  a  wedge  shock  wave  and  flow  field).  Also,  at 
low  angles  of  attack,  the  pressure  data  indicate  that  viscous  effects  have  become 
significant.  These  effects  tend  to  increase  the  value  of  the  similarity  parameter. 
Analytic  methods  for  determining  these  effects  on  delta-wing  flow  fields  are  not 
yet  available.  The  analysis  presented  here  neglects  viscous  effects. 

b.  Real -Gas  Effect  on  the  Similarity  Parameter 

The  value  of  the  similarity  parameter  on  delta  wings  is  increased  by  real- 
gas  effects.  This  increase  is  due  to  a  reduction  in  the  amount  of  mass  that 
leaves  the  lower  surface  over  the  leading  edges.  An  estimate  of  the  increase 
in  S  is  obtained  using  the  following  analysis. 


Assume  that  a  sharp  flat  delta  wing  at  angle  of  attack  a  generates  a  flat 
shock  wave  at  an  angle  0  to  the  free-stream  flow  as  illustrated  in  the  sketch 
below 


Section  A -A 


Also  assume  that  the  density  p2  in  the  shock  layer  is  constant  and  that  the 
lower  surface  velocity  u2  in  the  streamwise  direction  is  constant  and  of  a  magni¬ 
tude  equal  to  that  at  the  centerline.  Assume  that  the  flow  in  the  plane  normal  to 
the  lower  surface  at  the  leading  edge  is  sonic  and  directed  outward.  This  model 
is  reasonable  at  high  angles  of  attack  where  the  delta  wing  resembles  a  blunt 
body  in  the  plane  normal  to  the  wing  surface. 
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The  mass  rate  entering  into  the  shock  layer  is 


min  =  <u»  sln  6  ) 


bf 


(B2) 


(B3) 


2  cos  (0 -at) 

The  mass  rate  leaving  at  the  trailing  edge  is 


mte  =  °2  u2  b  f 


The  mass  rate  leaving  at  the  leading  edge  is 

"he  =  2  b  f  (P*a’/cos  A  "  P2  V  (B4) 

where  the  asterisk  indicates  sonic  conditions.  Balancing  the  mass  flow,  we  get 

“in  =  “te  +  “le  <B5a> 

P„  u tu  slnflj'COg(0_a)  =  f> 2  u2  bf  +  |bc  (p*a*/cos  A  -  P.,  u2>  (B5b) 

Now  make  the  assumption  that  P  *  as  p  2.  For  purposes  in  this  analysis, 
this  assumption  will  intensify  the  expected  differences. 

Divide  the  above  equation  by  p  2  03  be 

P„  u„ 


t  sin  8 


=  1  + 


a* 


P2  “2  c  cos(e-o)  u2  cos  A 

Because  c /f  =  tan  (0-a)  and  u2/u,„scos  0/cos(0-a),  then 
P2  tan(Q-a) _ 1 _ 


(B6) 


(B7) 


P-  ta"  8  "  1  +  secA 

U2 

The  left-hand  side  of  the  above  equation,  by  definition,  is  the  similarity  para¬ 
meter  so  that 


S  = 


1  +  a*  a 

—  sec  A 


(B8) 
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A  more  accurate  relationship  for  S  was  developed  by  correlating  sharp  delta 
wing  flow  field  results  obtained  using  the  method  of  reference  36.  This  corre¬ 
lation  resulted  in  the  expression 


S  = 


1 


1  +  0,6  a; 

u 


sec  A 


(B9) 


where  u  is  the  centerline  velocity. 

The  factor,  C.  6,  in  equation  (B9)  takes  into  account  the  changes  in  velocity, 
density  in  the  shock  layer, and  the  shock-wave  profile  in  the  spanwise  direction. 
Equation  (B9)  is  compared  to  the  results  of  reference  36  in  figure  B3.  The 
symbols  shown  in  this  figure  represent  solutions  for  sharp  delta  wings  at  var¬ 
ious  angles  of  attack,  angles  of  leading  edge  sw'eep,  Mach  numbers,  and  ratios 
of  specific  heats,  Y ,  Independence  of  Y  indicates  that  the  correlation  is  appli¬ 
cable  to  real-gas  flows.  Because  a*  can  be  significantly  smaller  in  a  real  gas 
that  it  is  in  an  ideal  gas  at  high  angles  of  attack,  the  real-gas  value  for  S  can 
be  substantially  larger  than  the  ideal -gas  value.  Results  shown  in  figure  B4 
indicate  that  changes  in  S  are  significant  at  high  angles  of  attack,  but  the  result¬ 
ing  change  in  the  shock-wave  angle  is  small. 


At  low  angles  of  attack,  values  for  S  obtained  using  equation  (B9)  are 
slightly  higher  than  those  proposed  in  figure  B2.  At  low  angles  of  attack,  how¬ 
ever,  equation  (B9)  Is  not  expected  to  be  valid  because  sec  A  is  no  longer  a 
characteristic  parameter  of  the  flow  field.  Also,  real-gas  effects  on  S  at  low 
angles  of  attack  are  small.  S  is  obtained  directly  from  figure  B2  at  low  and 
moderate  angles  of  attack. 


The  ratio  of  the  differences  in  S  that  can  be  expected  for  real-gas  flows 
and  flows  at  v  =  1. 4  has  been  correlated  in  terms  of  the  leading  edge  sweep 
angle,  A  ,  angle  of  attack,  a  ,  and  the  edge  enthalpy,  ie,  in  Btu/lbm.  This 
correlation  is  represented  by  the  relationship 


Sreal 

^ideal 


1  +  0.385 


/A\ 

/  a  \ 

(700) 

\  60*  / 

l°8l0  4e  "  2*7  I 

(BIO) 


which  is  compared  with  analytical  data  In  figure  B5. 


2.  SOLUTION  FOR  THERMODYNAMIC  GAS  PROPERTIES 


The  method  used  to  obtain  gas  properties  along  the  centerline  of  the  sharp 
delta  wing  surface  consists  of  two  parts.  First,  shock-wave  inclination  and  the 
corresponding  properties  downstream  of  the  shock  wave  are  determined  in  the 
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plane  normal  to  the  delta-wing  surface  at  the  centerline.  Second,  the  proper¬ 
ties  at  the  delta-wing  surface  are  calculated  by  using  the  similarity  parameter 
In  evaluating  the  surface  pressures.  Each  of  the  processes  is  discussed  in  sub¬ 
sequent  sections. 

a.  Shock-Wave  Solution 

The  solution  across  the  shock  wave  is  obtained  by  simultaneously  satisfying 
the  standard  shock  relationships 


P  u 

cn  o 
O 

P  +  O  U 


P  U 
2  2 


P  u 
2  2 


i  +  ^=i 

«  2  2  2 


2 


(Bll) 

(B12) 

(B13) 


and  the  expression  for  the  similarity  parameter  defined  by  equation  (Bl).  The 
symbols  p  ,  P,  and  i  represent  density,  pressure,  and  enthalpy,  respectively, 
u  is  the  velocity  normal  to  the  shock  wave,  and  the  subscripts  *  and  2  refer 
to  conditions  in  the  tree-stream  and  immediately  downstream  of  the  shock  wave, 
respectively. 

b.  Compression  Technique  Through  the  Shock  Laver 

After  the  shock-wave  angle  is  determined  and  the  gas  properties  downstream 
of  the  shock  wave  evaluated,  compression  through  the  shock  layer  to  the  surface 
of  the  delta  wing  must  be  computed.  For  conical  flows,  the  compression  takes 
place  isentropically.  To  evaluate  this  compression,  however,  one  of  the  gas 
properties  at  the  surface  must  be  known. 


Surface  properties  for  the  delta  wing  can  be  determined  analytically  using 
the  method  of  characteristics,  but  because  such  a  method  is  very  difficult,  an 
approximate  technique  was  developed. 

This  technique  is  based  on  a  correlation  between  the  pressure  ratio  (Pe/Ps) 
on  a  delta  wing,  the  pressure  ratio  on  a  cone,  and  the  similarity  parameter. 

The  similarity  parameter  S  was  developed  because  a  correlation  equation  such 
as  the  one  proposed,  which  shows  real-gas  and  leading-edge  sweep-angle  ef¬ 
fects  on  a  surface  gas  property,  was  desired.  A  correlation  equation  in  terms 
of  these  three  parameters  is  expressed  as 

(P  /P  )_  ..  =1+2  f  CP  /p  )  _  -  l]  (1  -  S)  (B14) 

e  s  Delta  wing  L  e  s  Cone  J  '  ' 
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Equation  (B14)  is  justified  using  the  following  physical  reasoning  and  analytic 
solutions. 

Figure  BG  shows  a  ratio  of  pressures  at  the  body  surface  and  at  the  shock 
wave  for  sharp  cones,  as  a  function  of  Mach  number  normal  to  the  shock  wave. 
The  curves  shown  on  this  figure  represent  solutions  using  the  method  of  charac¬ 
teristics.  Because  sharp  delta  wing  flows  at  high  angles  of  attack  are  also  con¬ 
ical,  it  can  be  expected  that  the  value  of  the  pressure  ratio  (P  /P  ),  for  the 
sharp  delta  wing  will  be  of  the  same  order  as  that  for  the  sharp  cone.  In  addi¬ 
tion,  pressure  ratios  obtained  using  the  method  of  reference  36  can  be  correlated 
in  terms  of  S  to  those  obtained  for  a  cone  using  the  method  of  characteristics. 

A  comparison  of  equation  (B14)  and  analytic  delta  wing  data  is  presented  in 
figure  B7  for  an  ideal  gas  with  >  =  1.4.  Values  for  (P  /P  )  in  a  real-gas 
flow  are  obtained  using  this  correlation  by  first  determining8  (P  /P  )„  and 
S  in  a  real-gas  flow.  e  s  Cone 

3.  PRESSURE  COMPARISON 

Surface  pressures  at  the  centerline  of  a  sharp  delta  wing  obtained  using  the 
method  described  above  are  compared  with  pressures  calculated  using  the  Boeing 
modified  Newtonian  method  in  figure  B8.  Comparisons  are  shown  for  a  70-de- 
gree  swept  sharp  delta  wing,  a  wedge,  and  a  sharp  cone  as  a  function  of  the 
surface  deflection  angle.  Differences  between  the  two  methods  result  in  pressure 
differences  on  the  delta  wing  on  the  order  of  2-4%  in  a  real  gas  and  7%  in  an 
ideal  gas. 

4.  STREAMLINE  DIVERGENCE  PARAMETER 

Sharp  delta  wing  centerline  heating  is  influenced  by  the  cross-flow  pressure 
gradient  at  the  centerline.  A  method  for  determining  this  influence  in  an  ideal 
gas  was  presented  in  reference  18  in  terms  of  a  streamline  divergence  parame¬ 
ter,  n,  defined  as 

-  /d  6  \ 

"  “\d0  i  (B15) 

0=0 


where  0  is  the  local  outflow  angle,  0  is  the  local  ray  angle  and  the  subscript, 
0-0,  denotes  that  evaluation  of  the  derivative  is  made  at  the  centerline.  The 
streamline  divergence  parameter  appears  explicitly  in  the  heat -transfer  rela¬ 
tionships  at  the  centerline,  as  was  shown  in  Appendix  A.  In  this  report,  results 
of  reference  18  have  been  extended  to  include  real-gas  effects. 

A  correlation  for  n  at  high  angles  of  attack  was  obtained  from  analytic 
solutions  of  refenence  36  and  is  shown  in  figure  B9.  This  correlation  can  be 
represented  by  the  expression 
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(BIG) 


n  =  0.53 


A 


] 


where  a*  is  the  sonic  velocity  at  the  leading  edge,  u  is  the  centerline  velocity, 
A  is  the  sweep  angle,  and  y  is  the  ratio  of  specific  heats.  The  parameter 
a*/u  sec  A  is  the  same  parameter  that  appears  in  the  correlation  for  S  at 
high  angles  of  attack. 

The  variation  in  n  is  illustrated  schematically  in  terms  of  angle  of  attack 
in  the  following  sketch. 


In  the  high  angle  of  attack  region,  equation  (BIG)  is  applicable.  Values  for  n 
in  this  region  are  shown  for  typical  wind-tunnel  and  flight  conditions  at  a  Mach 
number  on  the  order  of  20.  The  value  for  n  is  smaller  for  flight  because  the 
quantity  a*/u  in  equation  (BIG)  is  smaller  in  night.  At  low  angles  of  attack,  n 
is  calculated  by  the  ideal-gas  method  proposed  in  reference  18.  That  method 
consists  of  a  two-dimensional  flow  analysis  near  the  leading  edges  for  the  angle 
of  attack  region  where  the  shock  wave  is  attached  to  the  leading  edge.  Because 
shock  waves  are  attached  only  at  very  low  angles  of  attack,  the  method  is  strictly 
valid  only  in  that  region.  Typical  values  are  shown  for  this  region  in  the  left- 
hand  portion  of  the  preceding  sketch. 

In  the  intermediate  angle  of  attack  region,  no  theoretical  solution  for  n 
is  available.  In  reference  18,  an  interpolation  between  the  solutions  at  low  and 
high  angles  of  attack  was  proposed.  The  high  angle  of  attack  solution  is  that  of 
reference  36  for  an  ideal  gas  with  T  =1.4.  The  low  and  high  angle  of  attack 
solutions  were  "fitted"  to  one  another  so  that  the  interpolation  resulted  in  the 
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real-gas  effects  by  the  equation  shown  below.  correctea  lor 
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rection  is  made  The  _ .•  •  ucgrees  angle  of  attack,  no  cor- 

"weightine"  to  !he  II  function  ,n  equation  (BIT)  provides  a  satisfactory 
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Figure  Bl:  SIMILARITY  PARAMETER  FOR  SHARP  UNYAWED  CONES 


I 


Figure  B2:  SIMILARITY  PARAMETER  FOR  SHARP  UNYAWED 

DELTA  WINGS 
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Angles  of  attack  range  from  50°  to  70° 


I 
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Figure  B3:  CORRELATION  OF  SIMILARITY  PARAMETER  FOR  SHARP  DELTA 

WINGS  AT  HIGH  ANGLES  OF  ATTACK 


Delta  wing  surface 


0  ,  degrees 


~ldeal  gas  shock 

'-Real  gas  shock  with  uncorrected  S 
-Real  gas  shock  with  corrected  S  ,  eq.  (B1 
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Figure  B4:  EFFECT  OF  CHANGES  IN  SIMILARITY  PARAMETER  ON 
CENTERLINE  SHOCK  WAVE  GEOMETRY  FOR  SHARP 
DELTA  WINGS 
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Figure  B5:  REAL-GAS  CORRECTION  TO  SIMILARITY  PARAMETER 
FOR  SHARP  DELTA  WINGS 
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P  -  Surface  pressure  by  present  technique 

PmN  ■  Surface  pressure  by  Boeing  modified  Newtonian  method  (reference  2) 
Ideal  gas,  Moo  =  21.86,  Y  =  1 .4 

- R«°l  9«,  Voo  =  20,700  fps,  altitude  =  240,000  ft  (Moo=  21.86) 


Figure  B8:  COMPARISON  OF  SURFACE  PRESSURE  WITH 
MODIFIED  NEWTONIAN  VALUES 
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See  notes  on  figure  B3 
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Figure  B9:  CORRELATION  OF  DIVERGENCE  PARAMETER  FOR  SHARP 
DELTA  WINGS  AT  HIGH  ANGLES  OF  ATTACK 


=  0.53  (T) 


A,  deg 


raf.  18 


Figure  BIO:  DIVERGENCE  PARAMETER  FOR  SHARP  DELTA  WINGS 
AT  LOW  AND  INTERMEDIATE  ANGLES  OF  ATTACK 


APPENDIX  C 


METHOD  OF  SPANWISE  FLOW- FIELD  SOLUTION  FOR 
SHARP  UNYAWED  DELTA  WINGS 


Spanwise  flow-tield  distributions  are  developed  lor  the  low  and  the  high  angle 
of  attack  regions  for  an  infinitely  thin  sharp  unyawed  delta  wing  in  an  inviscid 
flow  field.  In  the  intermediate  angle  of  attack  region,  an  interpolation  between 
the  low  and  the  high  angle  of  attack  distributions  is  proposed  and  discussed. 

1.  LOW  ANGLE  OF  ATTACK  SOLUTIONS 

At  zero-degrees  angle  of  attack  in  an  inviscid  flow  field,  the  streamlines 
cross  the  leading  edges  of  a  thin  sharp  delta  wing  and  always  remain  parallel 
to  the  centerline  of  the  wing.  The  thermodynamic  gas  properties  and  the  stream¬ 
line  direction  on  the  lower  surface  are  everywhere  constant.  With  a  slight 
increase  in  the  angle  of  attack,  however,  the  flow  field  begins  to  undergo  a 
change.  At  the  leading  edge,  the  streamlines  are  deflected  slightly  away  from 
the  centerline  by  an  angle  0  *  as  illustrated  in  the  sketch  below: 

1 

0/P 


o 

When  the  shock  wave  is  attached  to  the  leading  edges,  0*  can  be  determined  from 
the  relationship 


where  Vn  and  Vt  are  the  velocity  components  on  the  surface  perpendicular  and 
parallel  to  the  leading  edge,  respectively.  The  velocity  component  Vn  is 
determined  using  the  oblique-shock  method  with  effective  values  for  the  Mach 
number  and  the  flow-deflection  angle.  The  effective  values  can  be  expressed 
in  terms  of  the  free-stream  Mach  number,  angle  of  attack,  a  ,  and  geome¬ 
tric  sweep  angle,  A  ,  by  the  equations  shown  below: 


6eff  =  tan'1  (tan  a  •  sec  A  )  (C3) 
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The  velocity  component  Vt  can  also  be  expressed  in  terms  of  Mm,  a  ,  and  A 
by  the  relationship 


vt  =  (M«  )(cos  <*  )(  sin  A) 


(C  4) 


At  the  centerline,  the  local  streamline  angle,  0  ,  is  zero  because  of  symme¬ 
try.  In  reference  18,  0  was  assumed  to  be  constant  within  a  region  near  the 
leading  edge.  This  region  was  bounded  by  the  outer  Mach  line  of  the  velocity 
vector  at  the  leading  edge  and  the  leading  edge  itself.  Between  the  inner  boundary 
of  this  region  and  the  centerline,  a  sinusoidal  variation  of  0  was  assumed.  In 
this  study,  0  was  assumed  to  vary  according  to  the  relationship 


JL 

6* 


(C5) 


where  £  is  on  the  order  of  unity  or  less.  A  value  C  =  0.5  approximates  the 
distribution  proposed  in  reference  18  reasonably  well.  It  was  shown  in  Appendix 
A  that  at  low  angles  of  attack  the  value  of  C  does  not  affect  the  heat-transfer 
results  appreciably. 

Pressure  on  the  surface  can  be  calculated  by  using  an  oblique-shock  method 
at  the  leading  edge  and  at  the  centerline  for  these  cases  where  the  shock  wave 
is  attached  to  the  leading  edge.  Typical  values  of  the  ratio  of  leading-edge  pres¬ 
sure  to  centerline  pressure  are  shown  in  the  following  sketch. 


The  spanwise  distribution  of  pressure  must  have  a  zero  slope  at  the  centerline 
because  of  symmetry.  At  the  leading  edge  the  slope  must  also  be  zero  because 
of  the  assumption  of  oblique-shock  flow  in  that  region. 
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2.  HIGH  ANGLE  OF  ATTACK  SOLUTIONS 

At  high  angles  of  attack,  spanwise  flow-field  solutions  can  be  obtained 
using  the  analytical  method  of  reference  36.  In  that  reference,  it  is  shown 
that  the  How  in  the  plane  perpendicular  to  the  delta  wing  surface  can  be  treated 
as  flow  over  a  blunt  body  having  a  cross-flow  velocity  at  the  leading  edge  equal 
to  the  local  sonic  velocity.  Because  of  conical  symmetry,  all  properties  along 
ray  lines  emanating  irom  the  apex  are  constant  and  the  calculation  involves  only 
two  dimensions.  All  streamlines  along  the  surface  appear  to  emanate  at  the 
apex  of  the  wing  as  shown  on  the  left  in  the  sketch  below.  The  corresponding 
distribution  of  the  local  streamline  angle  is  shown  on  the  right. 


In  the  subsequent  development  of  spanwise  distributions,  knowledge  or  the 
value  of  the  streamline  angle,  0*.  at  the  leading  edge  is  essential.  A  correla¬ 
tion  for  0*  from  the  analytic  solutions  has  been  made  in  terms  of  the  parameter 
a*/u  and  0,  where  a*  is  the  local  sonic  velocity  at  the  leading  edge,  u  is  the 
centerline  velocity,  and  0  is  defined  as  (90°  -  A).  This  correlation  is  shown 
in  figure  Cl  for  several  values  of  M^,  a  ,  A,  and  Y , 


0* 


tan 


0 


(C6) 


The  parameter  0*  serves  as  both  a  normalizing  and  a  correlating  factor  for  the 
local  pressure  and  the  local  streamline  angle.  The  spanwise  pressure  distribu¬ 
tion  is  required  to  determine  the  local  thermodynamic  gas  properties,  whereas 
the  local  streamline  angle  is  essenti.,1  in  determining  the  amount  of  the  local 
flow  divergence.  Distributions  Tor  these  parameters  are  shown  in  figures  C2 
and  C3  in  terms  of  the  parameter,  n/[e*  sec  A],  The  correlations  which  these 
distributions  are  based  on  are  shown  in  figures  C4  and  C5. 

The  distributions  for  n/ [o*  sec  A]>  1  are  shown  by  dashed  lines. 
Physically,  it  is  unlikely  that  values  of  n/[e*  sec  A  ]  >  1  can  occur.  This 
conclusion  may  be  reached  by  recognizing  that  the  local  slope  of  the  streamline 
angle  distribution  must  decrease  in  the  spanwise  direction  when  n/[0*  sec  A  ] 

>  1.  Correspondingly,  the  pressure  must  increase  in  the  spanwise  direction 
which  is  in  opposition  to  the  results  obtained  by  blunt-body  analogy.  Solutions 
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with  n/  [  0  *  sec  A  ]  >  1  indicate  that  some  of  the  assumptions  in  the  method 
are  not  satisfied.  However,  we  are  using  the  best  available  method  and  it  would 
be  a  major  task  to  improve  on  these  solutions. 

3.  INTERMEDIATE  ANGLE  OF  ATTACK  SOLUTIONS 

The  lower-surface  flow-field  results  at  high  and  low  angles  of  attack  have 
been  discussed  in  previous  sections  and  may  be  summarized  as  follows.  The 
streamline  angle  at  the  leading  edge  increases  from  the  value  zero  at  zero-de¬ 
grees  angle  of  attack  to  a  value  greater  than  0  a1,  high  angles  of  attack.  Physi¬ 
cally,  at  low  angles  of  attack,  the  ilow  enters  over  the  leading  edge  and  is  dis¬ 
charged  over  the  trailing  edge.  At  high  angles  of  attack,  the  surface  flow  behaves 
as  if  it  originated  at  the  apex  of  the  wing.  It  flows  downstream,  diverging  away 
from  the  centerline  and  eventually  discharges  over  the  leading  edge.  The  stream¬ 
line  angle  and  pressure  distributions  can  be  expressed  in  terms  of  the  local  ray 
angle,  <J> ,  and  at  high  angles  of  attack  are  uniquely  defined  by  an  additional 
parameter  n/[  0*  sec  A]. 

A  spanwise  distribution  of  the  local  streamline  angle,  0,  may  be  construc¬ 
ted  in  the  intermediate  angle  of  attack  region  using  the  preceding  development. 

By  knowing  both  the  value  and  the  slope  of  the  distribution  at  the  centerline  and 
the  value  at  the  leading  edge,  a  distribution  profile  between  the  leading  edge  and 
the  centerline  can  be  assumed.  Values  for  the  slope,  n,  and  the  streamline 
angle,  0*.  may  be  obtained  using  correlations  developed  previously  in  the  text. 
The  streamline  angle  at  the  centerline  is  zero  because  of  symmetry. 

At  some  intermediate  angle  of  attack,  0^,  the  streamline  angle,  0*,  at 
the  leading  edge  is  equal  to  the  leading-edge  ray  angle,  0.  For  all  angles  of 
attack  greater  than  the  angle  CK^,  Ilow  is  "expelled"  over  the  leading  edge  (i.e., 

0*  >  0  )•  For  all  G  <  flow  enters  the  lower  surface  by  crossing  the  leading 
edge  (i.e.,  0*  <  0).  The  streamline  pattern  and  the  corresponding  streamline 
angle  distribution  for  the  case  where  0*  -  0  are  illustrated  in  the  sketch  below. 


The  slope  n  at  the  centerline  is  less  than  unity. 

With  increasing  ang.e  of  attack,  outflow  occurs  over  the  leading  edges 
(i.e.,  0*  >  0).  The  correlation  for  n  indicates  that  the  value  of  n  remains 
below  unity  for  a  range  of  angles  of  attack  above  G1.  The  angle  of  attack  at 
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which  n  1  can  be  denoted  by  0(2-  In  the  angle  of  attack  region  between  0^ 
and  6*  >  f)  and  n  <  1.  The  streamline  pattern  and  the  streamline  angle 

distribution  in  the  region  ttj  <  Oc  <  Ctr,  are  typified  in  the  sketch  below. 


The  flow  field  has  a  "parting  line"  which  is  defined  by  the  ray  angle,  8  <j>  . 

The  flow  outboard  of  this  parting  line  spills  over  the  leading  edges  of  the  wing. 
The  flow  on  the  inboard  side  passes  over  the  trailing  edges.  A  particle  at  the 
location  where  8  =  <f>  remains  at  that  angular  location  with  increasing  distance 
from  the  apex.  In  order  for  the  streamlines  to  assume  a  pattern  as  shown  in 
the  sketch  above,  pressure  must  be  a  maximum  at  the  parting  line.  The  slope 
at  the  centerline  must  be  zero  and  the  value  at  the  leading  edge  must  be  on  the 
order  of  the  pressure  which  corresponds  to  the  local  sonic  pressure. 

The  change  in  location  of  the  "parting  line"  with  angle  of  attack  is  shown 
schematically  in  the  following  sketch: 


At  the  angle  of  attack  a^,  the  dividing  streamline  is  located  at  the  leading  edge, 
i.e.,  8*  /3  .  With  increasing  angle  of  attack,  this  dividing  streamline  shifts 

inward  and  eventually  reaches  the  centerline  at  an  angle  of  attack,  a0.  Above 
a2,  and  below  a^,  there  is  no  dividing  streamline.  Between  the  limits  otj 
and  tt„,  however,  the  dividing  streamline  plays  an  important  part  in  determining 
the  flow-field  and  heat-transfer  distributions. 

4.  FLOW-FIELD  REGIONS 

A  composite  of  all  types  of  flow-field  patterns  is  shown  schematically 
in  figure  C6  in  terms  of  angle  of  attack  and  Mach  number.  In  this  figure,  the 
first  flow-field  region  is  bounded  by  a  =  0  and  tne  angle  of  attack  at  which  the 
leading-edge  streamline  angle,  8*,  coincides  with  the  leading  edge  (i.e.  , 

0*  f}).  Within  this  region,  flow  passes  over  the  leading  edge  onto  the  lower 
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surface.  This  condition  is  illustrated  by  the  lower  three  sketches  on  the  right- 
hand  portion  of  this  figure.  Sketch  1  illustrates  the  lower -surface  streamline 
pattern  at  a  0,  sketch  2  shows  a  typical  pattern  in  this  flow-field  region, 
and  sketch  .3  indicates  the  pattern  when  0*  =  fi  ■  The  next  region  has  an  upper 
bound  defined  by  (d0/d01^,=o  =  1,  a  lower  bound  defined  by  0*  =0,  and  a  flow 
field  characterized  by  a  "parting  line"  or  dividing  streamline.  This  condition 
is  illustrated  in  sketch  4.  At  the  upper  boundary,  the  parting  line  has  reached 
the  centerline  of  the  lower  surface  and  all  flow  emanates  from  the  apex  of  the 
wing  and  away  from  the  centerline.  The  third  and  fourth  regions  occur  at  high 
angles  of  attack  and  are  characterized  by  the  condition  0  ><f>  for  all  spanwise 
locations  other  than  the  centerline.  The  third  region  corresponds  to  flows 
which  are  everywhere  supersonic  on  the  lower  surface.  It  is  in  this  region 
where  the  method  of  reference  36  is  valid.  The  fourth  region  occurs  at  very 
high  angles  of  attack  where  some  portion  of  the  lower -surface  flow-field  velo¬ 
city  is  sonic  or  subsonic.  The  lower  bound  of  this  region  is  defined  by  the 
condition  where  the  Mach  number  at  the  centerline  is  unity.  At  higher  angles 
of  attack,  the  total  velocity  component  becomes  subsonic  at  the  centerline  and 
moves  the  sonic  line  further  outboard.  At  90 -degrees  angle  of  attack,  the  sonic 
line  would  coincide  with  the  leading  edge.  The  subsonic  region  of  flow  is  out¬ 
side  the  limits  of  the  present  analysis. 

Streamline  angle  profiles  over  the  lower  surface  of  the  sharp  della  wing 
are  shown  schematically  in  figure  C7.  These  profiles  have  been  drawn  to 
correspond  to  the  flow-field  regions  shown  on  figure  C6.  At  a  =  0,  the  stream¬ 
line  angle,  0,  is  zero,  and  all  streamlines  are  parallel  to  the  centerline.  This 
is  depicted  in  sketch  1.  As  the  angle  of  attack  is  increased,  the  streamlines 
are  deflected  away  from  the  centerline  as  they  cross  the  leading  edge.  This  is 
shown  by  the  streamline  angle  profiles  pertaining  to  sketch  2.  The  situation 
where  the  leading -edge  streamline  angle,  0*,  is  equal  to  the  leading -edge 
sweep  angle,  /3,  is  depicted  by  the  distribution  corresponding  to  sketch  3.  At 
this  point,  the  dividing  streamline  coincides  with  the  leading  edge.  As  the 
angle  of  attack  of  the  delta  wing  is  increased  0*  becomes  larger  than  fi  and  the 
lower  surface  assumes  a  streamline  angle  distribution  like  that  shown  for 
sketch  4.  The  dividing  streamline  has  shifted  away  from  the  leading  edge  and 
is  located  at  9)g=0.  With  increasing  angle  of  attack,  the  dividing  line  shifts 
more  and  more  toward  the  centerline.  The  point  at  which  the  dividing  stream¬ 
line  has  just  reached  the  centerline  is  defined  by  (d0/d<^  0  =  1.  The  flow 
field  for  this  situation  is  typified  in  sketch  5.  At  still  higher  angles  of  attack, 
(d0/d0)0_Q  >1  and  the  streamline  angle  distribution  is  like  that  shown  for 
sketch  G. 

Spanwise  pressure  profiles  corresponding  to  the  flow  regions  discussed  in 
the  previous  paragraphs  are  shown  in  figure  C8.  At  zero  degrees  angle  of 
attack,  no  flow  divergence  is  encountered  and  the  pressure  distribution  is 
constant.  At  low  and  intermediate  angles  of  attack  the  pressure  increases  in 
the  spanwise  direction.  This  is  illustrated  by  the  distributions  pertaining  to 
sketches  2  and  3  in  the  figure.  For  a  flow  field  with  a  parting  streamline. 
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pressure  is  at  a  maximum  at  the  parting  streamline.  This  is  indicated  by  the 
distribution  pertaining  to  sketch  4.  For  this  case,  the  pressure  at  the  leading 
edge  is  on  the  order  of  the  local  sonic  pressure.  At  high  angles  of  attack,  the 
pressure  distributions  are  typified  by  those  pertaining  to  sketches  5  and  6  in 
the  figure.  The  pressures  at  the  leading  edge  are  equal  to  the  local  sonic 
pressures. 

The  location  of  the  dividing  streamline,  defined  by  the  ray  angle  8  =<f> , 
is  sensitive  to  the  angle  of  sweep,  the  free-stream  Mach  number,  and  the  total 
enthalpy.  This  is  depicted  in  figure  C9  for  a  delta  wing  of  70 -degrees  sweep 
at  free-stream  Mach  numbers  of  10  and  20.  The  important  feature  to  be  recog¬ 
nized  in  this  figure  is  the  marked  difference  in  location  of  the  dividing  stream¬ 
line  between  flows  of  different  total  enthalpies. 
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Figure  Cl:  CORRELATION  FOR  LEADING  EDGE  STREAMLINE  ANGLES 

AT  HIGH  ANGLES  OF  ATTACK 
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SPANWISE  STREAMLINE  ANGLE  DISTRIBUTION  AT 
HIGH  ANGLES  OF  ATTACK 


Figure  C2 
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Figure  C3:  SPANWISE  PRESSURE  DISTRIBUTION  AT 
HIGH  ANGLES  OF  ATTACK 
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n/[0*  sec  A] 

Figure  C4:  CORRELATION  FOP  SPANWISE  STREAMLINE  ANGLES 

AT  HIGH  ANGLES  OF  ATTACK 
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Figure  C5:  CORRELATION  FOR  SPANWISE  PRESSURE 
AT  HIGH  ANGLES  OF  ATTACK 
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Figure  C6:  LOWER  SURFACE  FLOW  FIELD  REGIMES 
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Figure  C7:  LOWER  SURFACE  STREAMLINE  ANGLE  PROFILES 
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Figure  C8:  SPANWISE  PRESSURE  PROFILES 
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